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Abstract 



■ The usual sphaleron bound and the statement of the impossibility of baryon production at a 
CD . second order phase transition or analytic cross-over are reformulated in the first part of the 

I paper as requirements of the expansion rate of the Universe at the electroweak scale. With 

Q>^ ' an (exact or effective) additional contribution to the energy density scaling as 1/a^, which 

■ dominates until just before nucleosynthesis, the observed baryon asymmetry may be produced 
at the electroweak scale in simple extensions of the Minimal Standard Model, even in the case 

^ that the phase transition is not first order. We focus our attention on one such cosmology, in 

Q_i' which the Universe goes through a period termed kination in which its energy is dominated by 

^ • the kinetic energy of a scalar field. The required kinetic energy dominated modes can occur 

D ! either as a field rolls down an exponential (or steeper) potential, or in the oscillation of a field 

' about the minimum of a steep power-law potential. We implement in detail the former case with 

a single exponential field first driving inflation, and then rolling into a kinetic energy dominated 
^ . mode. Reheating is achieved using an alternative to the usual mechanism due to Spokoiny Q, 

^ I in which the Universe is 'reheated' by particle creation in the expanding background. Density 

perturbations of the magnitude required for structure formation may also be generated. We 
show that the analogous model for the power-law potential cannot be consistently implemented. 
In models with inflation driven by a second field and the usual mechanism of reheating (by 
decay of the inflaton) the required kinetic energy dominated cosmology is viable in both types 
of potential. 
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1 Introduction 



Nucleosynthesis provides a role model for electroweak baryogenesis to whose impressive heights it can still 
only aspire. The great attraction of the idea that the baryon asymmetry of the Universe (BAU) may have 
been created at the electroweak epoch lies in the possibility that one day an ab initio calculation to rival 
that of nucleosynthesis may be possible, and that it will give a definitively positive or negative answer. 
Rather than simply providing an alternative to scenarios for baryon creation at the GUT scale, it has the 
fundamental interest of relying on physics at a scale directly accessible to experiments. We can realistically 
hope to know the correct theory of physics at the electroweak scale, in particular the structure of the CP 
violating and symmetry breaking sectors. Just as in nucleosynthesis it is then a question of putting this 
theory in an expanding universe and finding the output. Electroweak baryogenesis however faces more 
substantial obstacles on the road to a reliable calculation than did nucleosynthesis, e.g. the determination 
of the baryon asymmetry involves all the details of departure from equilibrium at the phase transition 
(if there is one), the crucial baryon number violating processes arising from the chiral anomaly at finite 
temperature involve many difficult and still unresolved questions, etc. Much progress has however been 
made, recently in particular using lattice methods to study the phase transition ||2| , and the problems seem 
not to be insurmountable. 

The approach of this paper is somewhat orthogonal to the direction of investigation of most work on 
electroweak baryogenesis. Rather than investigating some aspect of the particle physics, we consider the 
cosmological side of the problem. The standard and indeed most natural assumption about cosmology 
at the electroweak epoch is that it is what one gets by the simplest backward extrapolation from nucle- 
osynthesis: a homogeneous and isotropic radiation dominated universe. In nucleosynthesis the assumption 
of such a universe is relaxed to place limits on, for example, the contribution of a magnetic field or of a 
cosmological "constant" to the energy density. In this paper we ask the analogous question of electroweak 
baryogenesis: how is the standard scenario for production of the observed baryon asymmetry at the elec- 
troweak epoch affected if we consider cosmologies other than the standard one? And are there simple 
alternative cosmologies which lead to significantly different results for electroweak baryogenesis? 

The same sort of question has been previously addressed in the context of calculations of the relic 
densities of weakly interacting particles in work of Barrow |^ and Kamionkowski and Turner (KT) 
These relic densities depend on the temperature at freeze-out which occurs (approximately) when the 
annihilation cross-section of the particular species drops below the expansion rate of the Universe. Barrow 
discussed the particular case of a non-anisotropic universe, in which the average (volume) expansion rate 
which determines the freeze-out has an extra component driving it which scales as (a is the average 
scale factor). Consistent with the requirement of radiation domination at nucleosynthesis, the expansion 
rate can thus be very much greater in the anisotropic universe prior to nucleosynthesis when dark matter 
relics typically freeze-out (T ~ lOOMey), and the requirement that such a particle be the cosmological 
dark matter may in principle place a bound on the anisotropy. The important idea - that relic densities 
can provide a probe of the Universe prior to nucleosynthesis, which might be other than the standard 
radiation dominated one - was considered in a more general way by KT, who discussed the case of an 
anisotropic universe, as well as various others, including a Brans-Dicke-Jordan theory of gravity. In this 
latter case the eff'ect can also be modelled as an extra contribution to the energy density scaling as 
producing in the same way a speeded up expansion rate before nucleosynthesis without violating the 
nucleosynthesis constraints. KT also mention an example (which they describe as "exotic") of a scalar 
field (j) which oscillates in the minimum of a potential , for which the energy density scales as a~^"/ ("-+2) ^ 
i.e. faster than radiation for n > 4. Again the energy in such a mode can contribute significantly before 
nucleosynthesis without disrupting the latter. As discussed in the relevant feature of this model is that 
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it is the kinetic energy of the scalar field which gives the dominant contribution to the energy density of 
the Universe. As well as the oscillating mode of the power-law potential, the scaling applies to a scalar 
field rolling down a simple exponential potential. Rather than being exotic (compared to the models which 
^] and Q focus on), such models are minimal in the sense that they leave Einstein gravity intact and are 
consistent with the inflationary explanation of the homogeneity and isotropy of the Universe. In this paper 
we construct and study in detail a model for each of the two cases in which the single field (exponential or 
power-law) both inflates the Universe and then rolls into the kinetic energy dominated mode. Reheating is 
achieved using a simple alternative model of reheating proposed by Spokoiny several years ago In the 
power-law potential density perturbations are produced which are too large and the model is not viable. 
Furthermore the coherent oscillating mode is unstable to decay due to parametric resonance. We also 
discuss less constrained viable models in which the infiaton is a different field and reheating proceeds in 
the usual way (by decay of the condensate). The required potentials do in fact arise in many particle 
physics models: Power-law potentials have been discussed, for example, in the context of supersymmetry 
motivated inflationary models in |^, 0]. (The lower order terms can be excluded by imposing a discrete 
Zn/2+i symmetry on the superpotential). Exponential potentials arise quite generically in theories involving 
compactified dimensions, such as supergravity and higher dimensional theories of gravity (for specific 
examples, see [^, ^, |l^). The latter are also interesting in that they can play an important (potentially 
observable) role in the late-time cosmology of structure formation |11, 

In the first part of the paper we address the question of how the expansion rate of the Universe affects 
the baryon asymmetry produced at the electroweak scale, without reference to any particular cosmological 
model. The two distinct cases - a first order phase transition, and a second-order phase transition or a 
cross-over - are treated separately. In the first case the expansion rate enters (i) in determining when the 
transition occurs, since this depends on the cooling rate of the Universe below the critical temperature, 
and (ii) in determining the depletion of the baryon asymmetry produced by sphaleron processes (and hence 
in determining the sphaleron bound). In this case the baryon asymmetry is actually (at least in certain 
extensions of the standard model) produced on or near the bubble walls as they propagate through the 
plasma, and does not depend directly on the expansion rate. In the second case the expansion rate is the 
sole parameter which controls the departure from equilibrium and the baryon asymmetry which can be 
generated depends directly on it. 

In the second part of the paper we turn to the discussion of alternative cosmologies, first reviewing 
those considered in previous work and then turning to the detailed consideration of cosmologies dominated 
by the kinetic energy of a scalar field for a period between inflation and nucleosynthesis, concentrating 
in particular on the case when this phase (which, following ^ we term kination) persists until after the 
electroweak scale. In this case the expansion rate at the electroweak scale is increased, producing the 
modifications to calculations for electroweak baryogenesis discussed in the first part. As outlined above, 
the two types of models we consider are exponential potentials and power-law 0" (n > 4) potentials. For 
both cases we first discuss a one field model, in which the field both infiates the Universe and causes 
kination. Instead of decaying as in the standard explanation of reheating, the infiaton rolls into a kinetic 
mode and simply red-shifts away. The radiation created by the superluminal expansion of the Universe 
at the transition between the two phases thermalizes and comes to dominate the energy density of the 
Universe at a later time determined by the expansion rate at the end of inflation We show that in the 
exponential potential one can have (i) a transition to radiation domination as late as nucleosynthesis and (ii) 
thermalization of the radiation well prior to the electroweak scale, and further, (iii) density perturbations 
of the right amplitude for structure formation. In the power-law potential, however, we find that the 
requirement that the Universe become radiation dominated before nucleosynthesis leads to the production 
of density perturbations which are much larger than is consistent with observations. In any case the 
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oscillating mode in this model typically decays non-perturbatively (through parametric resonance) unless 
the self-couplings of the field are extremely small, and the energy does not stay in the kinetic mode for 
long enough. We conclude section 5 with a brief discussion of two fields models in which one field is the 
inflaton and reheats the Universe in the standard way, and the second field is the 'kinaton' which comes 
to dominate for a phase subsequent to inflation. In this case the second field can be either exponential or 
power-law (provided the couplings are such that decay by parametric resonance does not occur until after 
nucleosynthesis) . In the last section we discuss the implications of our results for the testability of theories 
of electroweak baryogenesis, and consider briefly other ways in which pre- nucleosynthesis cosmology might 
be probed. 

2 Dependence on the Expansion Rate 

In electroweak cosmology the assumption is generally made that the Universe is flat, homogeneous, 
isotropic, and radiation dominated. Hence all cosmological information is encoded in the expansion rate 
i^rad) which is given as a function of the plasma temperature T: 

H,,^ = , h=( ^1 ' (1) 

where g^: ~ 10^ is the number of relativistic degrees of freedom in the plasma and Mp = (SvrG)^^/^ ~ 
2.4 X lO^^GeV is the reduced Planck mass. The clean separation between the purely particle physics and 
cosmological calculations occurs because of the adiabaticity of the expansion. It is adiabatic because the 
timescale characterizing the expansion - Texpansion = 

Mp/{hT^) ~ 10^'^(100GeV/r)T-\ taking ~ 10^ - 
is much greater than the timescales associated with the thermalization processes which have typical rates 
~ a^T (where a > 1/50 for all the interactions well above lOOGeV). The phase transition can thus be 
studied using equilibrium methods - the expansion of the Universe enters only in determining the cooling 
rate and, hence, when the transition occurs (if it does). In general we could of course consider any cosmology 
at the electroweak scale, with the sole requirement that it be consistent with nucleosynthesis. We have 
no probe of the electroweak scale except that provided by electroweak physics, and methodologically it 
makes sense to ask how changing the standard assumption about the Universe at this phase affects the 
predictions of the remnants which result. Here we limit ourselves to relaxing only the assumption that the 
expansion rate is related to the temperature by (|l]). Instead we take 

H = H,^( — ) (2) 

where p is a number and the subscript 'ew' means the quantities are evaluated at some temperature char- 
acteristic of the electroweak phase transition. Using H = a/a this corresponds to the time dependence 
a oc t^^^ for the scale factor a. All our results concerning baryogenesis are, we will see, essentially inde- 
pendent of p because they depend only on temperatures very close to Tew We will treat Hew as a free 
parameter, only taking it to be such that the assumption of adiabaticity is valid, which allows it to be 
different form the standard value by orders of magnitude. We will review in section 5 some of the non- 
standard cosmologies which can be described by these assumptions. The model which we will discuss in 
detail is a homogeneous and isotropic universe dominated by a kinetic mode of a scalar field rather than 
by radiation. 

Baryogenesis, the creation of baryons from an initial zero baryon state, requires a departure from 
thermal equilibrium. In the Big Bang Universe this is provided by the expansion which causes the Universe 
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to cool. At the electroweak scale this cooling can lead to two very different effects, depending on whether 
the electroweak phase transition is of first order or not. Recall that at a first order phase transition as the 
Universe cools it becomes thermodynamically favourable for the system to be in the "broken" state. The 
"broken" and "unbroken" phases are separated by a potential barrier which decreases as the Universe cools. 
Once the barrier is low enough, the transition proceeds by the nucleation and propagation of true vacuum 
bubbles. This departure from equilibrium is characterized by timescales which are much shorter than 
that associated with the expansion. Almost all proposed mechanisms for baryogenesis at the electroweak 
scale make use of this dramatic departure from equilibrium, using the interaction between the plasma and 
propagating walls to generate the baryon asymmetry. In the case that the transition is second order or 
cross-over there is no such effect. Everything evolves continuously and the departure from equilibrium is 
controlled directly by the expansion rate. This usually leads one to conclude that anything but a first order 
phase transition is inimical to baryogenesis at the electroweak scale. Once one relaxes the assumption that 
the expansion rate is its standard radiation dominated value, this conclusion does not follow and needs to 
be examined more carefully. We will thus treat these two cases separately in some detail. 



3 First Order Phase Transition 

In this case the baryons are created as the bubbles of the true vacuum propagate through the false vacuum. 
The net effect of the propagation of the bubble through the medium in all proposed mechanisms is the 
creation of a flux of baryons into the broken phase. The expansion rate enters only indirectly through other 
parameters involved in this calculation - through the temperature at which the transition occurs which it 
determines, and, in certain regimes, through the bubble wall velocity. On the other hand, it enters directly 
in the determination of the amount of the created asymmetry which survives once it is in the broken phase. 
We consider these two dependences separately. 



3.1 Bubble Nucleation 

In this section we investigate how the bubble nucleation temperature depends on the expansion rate of 
the Universe. We also briefly discuss how the bubble wall velocity may depend on this parameter. As the 
Universe supercools below the critical temperature Tc, the fraction / filled by nucleated bubbles at a time 
t is given by (see Q and references therein): 

fit) = 1 - e-^W , A(t) = /* dt'^v^t - t'fn{t') , n = loT^e-'^^/^ , (3) 

where S'b is the bounce action, TZ the nucleation rate per unit volume, /q is a prefactor which is a slowly 
varying function of temperature of order one, given in more detail below, and v is the bubble wall velocity. 
Changing the integration variable to x = (Tc — T) /Tc and using the time-temperature relation t oc T^'p 
which follows from (|2[), one finds 

A(x) = Jo -f / dx'n - x'r P — r ; CXp °^ (4) 

3 ^\hJ Jo ^ ' \p{l-xy p{l-x')p) ^\ Tc{l-x')J 

where He is the expansion rate at Tc. We will see below that the nucleation temperature T^uci defined by 
A(xnuci) = 1 is always very close to the critical temperature so that we can take < x' < x <C 1 and 
expand to linear order in (^ to get 
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Keeping the first term in a derivative expansion of the bounce action about x, i.e. taking S'b(x') = Sh{x) + 
{dS]:,{x)/dx){x' — x) + 0{x' — x)"^, the integral can be performed with the assumption that xd{S\y/T) / dx{x 
Xc) ^ 1, and gives the nucleation temperature imphcitly as 
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where -ffrad — 1-2 x 10~^^(Tc/100GeV)rc is the expansion rate at Tc in the standard radiation dominated 
cosmology. 

To check the consistency of our assumptions and evaluate this expression to give the nucleation tem- 
perature, one must calculate the bounce action near the critical temperature in the particular model of 
interest. We consider the Minimal Standard Model (MSM) in the regime where it is described by the 
effective potential 
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with the one-loop ring improved values lO, |T4| 
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where vq = 246GeV, as = (47r)2e"2TB ~ 50, = {Tif'e''^'^'^ ~ 3.1, and 7e is Euler's constant. This 
treatment of the MSM is reasonably accurate up to mn ~ 60GeV, when nonperturbative effects become 
important. With this effective potential the critical temperature Tc is given by 
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and the latent heat L and surface tension a by 
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where defined by degeneracy of the minima T) = 1/(0, T)), is 

_2a- [4a2 _ 18At7(1 - (To/T)^)]^/' 



(10) 
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(12) 



with the boundary conditions ^{r = 0) = (j)^, d4>/dr{r = 0) = 0, and (j){r = oo) = (r = \x\ is the 
radial coordinate). Rather than solving this exactly (which is numerically expensive), or in the thin wall 



The (spherical) bounce action is given by 

S'b = 47r / r'^dr 
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approximation (which is inaccurate for strong phase transitions when cp^^ ~ T) , we will use an approximation 
for S'b developed in [p^ ]: 
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(13) 



which is valid for < A < 0.95. 

The prefactor Iq in (0) can be written as 
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where Kdyn is the dynamical prefactor as given in |15|. i? is the radius of the nucleating bubble, which 
can be estimated in the thin wall approximation to be i? ~ -Rnuci ~ 2a/Lxnuci, Xnuci = 1 — Tj^^ci/Tc, 
p = 7r^T^(7=K/30 is the energy density of the plasma with relativistic degrees of freedom. The one loop 
fluctuation determinant consists of the 'negative' mode A_ ~ O.O5(7^/^0(T), and we take /Cbubbie = 1 (for a 
more accurate value, see [16, |l7|] ). 

We have solved (^ numerically to find the nucleation temperature Tnucb using these values and ap- 
proximations. We also used mt = 175GeV, m^y = 81GeV, mz = 91GeV, and took the bubble wall velocity 
V = 0.4 [|l4| 1^. In Figure 1(a) we show a plot of Cnuci = (?c — T^vic\)/{Tc — Tq) against the logarithm of 
the expansion rate at the critical temperature, for a range of Higgs masses rriH- We see clearly that 
the usual result in a radiation dominated universe is qualitatively unchanged by varying the expansion 
rate over orders of magnitude: C,c is small, so the nucleation temperature is very close to the critical tem- 
perature, and much closer to than Tq. There is a small quantitative change, Cc varying by about 40% 
as the expansion rate varies over five orders of magnitude, but the change in absolute terms is tiny since 
ATnuci = -(Tc - ro)ACc < 5 X 10-3 

^nuci- Typically there is about the same change to C,c over this range 
at a fixed Higgs mass as is brought about by decreasing the Higgs mass by about 25GeV in a radiation 
dominated universe^ 

The results are easy to understand both quantitatively and qualitatively. Varying the bounce action 
( p!^ ) and assuming 1 — ^ < 0.2 (which is satisfied for most of the parameter space in figure 1(a)) so that the 
last term in the expression dominates, one obtains 5[S\JT) ~ — 2((5Cnuci/Cnuci)('S'b/r)0> and hence (using 
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Cnucl "Sb/T 

The bounce action 5b approximately halves, going from lOOT to 50T as the expansion rate changes from 
-f^rad ~^ lO^iJi-ad- Taking an average value for it in (|l5|) gives good agreement with the estimates we made 
above from the figure. Qualitatively the reason the expansion rate changes the nucleation temperature so 
little is that, as long as the Universe expands on a time-scale much longer than T^^^j, the transition always 

^This is a friction limited upper bound. The results here are of course not very sensitive to the details of the prefactor 
in the nucleation rate. As we will discuss below, this assumes that that the Universe is not reheated significantly by latent 
heat release, which is a reasonable approximation in the MSM. In the case that significant reheating occurs, so that bubble 
nucleation stops, the correct value of v in the early stages of nucleation would be the speed of sound Vs ~ l/VS at which the 
shock fronts propagate. 

^ Cnucl increases as the Higgs mass rnn decreases because the phase transition gets stronger, and therefore more supercooling 



^This estimate is just that obtained in the thin wall approximation in which Sb,thin oc 
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Figure 1 (a), ^^uci expansion rate 
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proceeds when the nucleation rate is very suppressed, where the bounce action is an extremely sensitive 
function of temperature. The nucleation temperature decreases as the expansion rate increases because 
the Universe must supercool more to attain a less suppressed nucleation rate. 

Such a small change to the nucleation temperature leads to minor changes to the quantities which 
determine the baryon asymmetry generated. We will see in section 3.2 that between Tc and Tq the VEV 
of the Higgs field changes by 50%, and its derivative with respect to T by about a factor of three. From 
figure 1(a) this would mean an increase in the VEV at nucleation of 1% or a little more per order of 
magnitude increase in the expansion rate. We would expect that this result will hold true in any typical 
electroweak model and not just the MSM in the regime we have studied it here. These minor changes to 
the VEV(s) and the other macroscopic parameters which determine the baryon asymmetry (bubble wall 
thickness, profile etc.) are essentially negligible in their effect on the baryon asymmetry generated. 

One condition must be attached to this conclusion, however: Other macroscopic effects can come into 
play as the bubbles propagate. When the propagating bubbles begin "bathing" in the hydrodynamic shock 
waves of the neighboring bubbles, i.e. when A(i) ~ 1 in (^), the plasma can heat up and slow down the 
propagation of the bubbles |18]. To determine how big this effect can be one compares the latent heat 
release L with the difference in the thermal energy density A/9 = 4:p{Tc — T^vlc\)/Tc between the nucleated 
phase and the unbroken phase. If L/ lS.p > 1 the system can reheat all the way back up to Tc. If such 
reheating occurs the main effect on baryon generation at the bubble walls is through the slowing down of 
the bubble walls. 

In figure 1(b) we show a plot of this ratio L/ lS.p as a function of the expansion rate in the MSM with the 
same values and range of Higgs mass as above. As the expansion rate increases the amount of re-heating 
decreases - simply a result of the increased supercooling. In all the parameter space the ratio is less than 
one, but of order one, so the effect of reheating may be significant, leading potentially to a slow-down of the 
bubble walls relative to the friction dominated regime. To determine the effect precisely is a very involved 
problem, and we will not attempt to tackle here the even more involved one of looking at the effect on 
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Figure 1(b). L/Ap vs the expansion rate 
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such slow-down of changing the expansion rate. We hmit ourselves to the qualitative observation that the 
amount of reheating decreases slightly as the expansion rate increases, and that modifying the expansion 
rate by orders of magnitude is likely to lead to only small corrections. 

The one-loop effective action which we use here becomes an increasingly poor approximation for niH > 
60GeV, and the latent heat is one of the quantities it estimates very inaccurately. Lattice studies Q have 
shown that there is a regime in the MSM where L/Ap > 1. In other models also this is certainly a condition 
we can envisage being satisfied. As mentioned above, in this case the bubble walls slow down to a final 
velocity determined by the expansion rate of the Universe, which has been estimated to be u ~ 
I p^ , ^ in a radiation dominated universe. These estimates show that this velocity depends on the expansion 
rate only through the combination HR where R is the average bubble radius, which is the radius at which 
5f/5R of (ID peaks, i.e. (fA/dR^ 0. A little algebra gives 



HR ~ 



3v. 



\dS/dT\ 



(16) 



where Vg = l/\/3 is the speed of sound. Now using the expression (^) above we see that \dS]^/dT\ oc 5b/T 
decreases by about 10% - and hence the final wall velocity increases by the same amount - per order of 
magnitude change in the expansion rate. If L/Ap > 1 for different expansion rates, this will be the only 
change to the final bubble velocity. Given the behaviour we have observed of L/Ap it is clear that in certain 
models increasing the expansion rate considerably relative to its radiation dominated value could make the 
difference between this final wall velocity and a (typically much larger) friction limited one corresponding 
to L/Ap < 1. In this particular case the change to the bubble wall velocity may not be so small. 

The determination of how the baryon asymmetry would be affected by such changes in the velocity of 
the bubble wall is a model dependent problem. The velocity dependence of the baryon asymmetry depends 
on what precise mechanism is operative, which depends on both the microscopic and macroscopic physics. 
At low velocities {v < 10~^) the result always goes to zero at least linearly, and at larger velocities the 
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most sensitive dependence is ~ Xjv^ . Using this dependence and assuming the greatest possible effect due 
to a change in the expansion rate (from an upper bound in the friction dominated regime v ~ 0.4 to the 
lower bound of the adiabatic 'complete reheat' regime with v ~ 10-3) would give a change in the calculated 
BAU by (at most) 10^ As we have discussed however, in most models the change will be much smaller 
and probably very small. A more detailed investigation of this question would however be required to draw 
stronger quantitative conclusions. 



3.2 Washout and the Sphaleron Bound 

The baryons created at the bubble walls are subject to decay after they enter the broken phase if the 
baryon number violating processes are not sufficiently suppressed. The requirement that this attenuation 
not reduce the created asymmetry to less than that required for nucleosynthesis leads to the sphaleron 



bound 1 19] in a radiation dominated universe. Here this bound will be restated as a requirement of the 
expansion rate of the Universe in a given theory. In the course of this discussion we will also draw attention 
to imprecisions in commonplace statements of the sphaleron bound (with the usual assumption of radiation 
domination) which can be of considerable importance. 

Since the time-scale associated with the baryon number violating sphaleron processes is much longer 
than the time-scale for thermalization processes, the baryon number after the completion of the electroweak 
phase transition is given by its local equilibrium value 

'Zi ^ -I 

where \xb and ha are chemical potentials for baryon number and the other charges Qa conserved on the 
relevant time-scale, i.e. exactly conserved, or violated at a rate slower than the baryon number violating 
processes. On the time scale over which the violation of baryon number is relevant, the system relaxes to 
equilibrium at a rate B = —Tspii{AF/T)AB, where Tsph is the rate per unit volume of sphaleron processes 
in which the Chern-Simons number A'cs changes by one unit, and AF is the free energy change per process. 
Since AB = NpAN^s per process {Np is the number of fermion families), we get 

B = -Nf-^„B (18) 
where we assume that other charges are defined so that they are conserved in these processes {AQa = 0). 



The sphaleron rate is given by [20| 



"■sph ^ ^ / 47r \ 2 _ ^ (j_ 



rsph = C5^e-— , ^sph = S — C = - -7;^MrA4otVrot/Csph (19) 

T'^ \awJ 2TTg(f){T) 

where a^j = g'^ /An ^ 1/29, B is a monotonically increasing function oi X/g^ = m^j/^m^ ranging between 
1.5 and 2.7 as \/g^ varies from to oo [^] and C is a temperature independent 'constant', given fully 
below (with uj_ ~ g(j){T) the frequency of the negative mode of the sphaleron, /Cgph the one loop fluctuation 
determinant, Vrot = Stt'^ a group volume factor, and Mtr and A^rot the number of translational and rotational 
degrees of freedom). 

The conserved charges Qa are just the primordial values of the exactly (or, in some cases, approx- 
imately) conserved charges in the electroweak model with which we are calculating. In scenarios for 
electroweak baryogenesis these are always taken to be zero. From ( p!7| ) it then follows that fis can be 
expressed in terms of B, so that (p^) becomes simply 



B = -OnTsphS , Tsph = (20) 
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where a„ is a number of order one whose precise value depends on the model and its corresponding set 
of charges Qa- In section 4 below we carry out the constraint calculation explicitly and find a„ ~ 0.4 for 
typical electroweak models. Integrating ( pO[ ) gives the baryon asymmetry -Breeze which survives to partake 
of nucleosynthesis: 



^freeze = -B(7b) CXp 



dtanTsphit) 



-B(7b)exp 



_ 1 f^-^^a^n^ 

^ Jo Vrj., 



(21) 



where B{T^) is the baryon asymmetry at the completion of the transition, at temperature Tb and H^, is 
the expansion rate at that time. As discussed in the previous section the appropriate value of Tb depends 
on the details of the of the phase transition and lies in the range £ [Tc, Tnuci]- To obtain the latter form 
of ( pl| ) we have used the time-temperature relation t oc which follows from (^). Changing variables to 
y = Ty^/T we can write the depletion factor D as 
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In- 



B 



freeze 



^ xQanCNlg 



(22) 



Over the range of integration the factor -Egph/^b in the exponential increases from its minimum value at 
y = 1, which is quite a large number ~ 30. This means that the dominant contribution to this integral 
comes from temperatures very close to Tb with y ~ 1 + 1/30. In fact we will see below that the rate of change 
of the VEV is typically large enough to narrow the range of temperatures which dominate the integral even 
more than this. The p dependence in the integral is therefore very weak and the only significant effect of 
the change in the expansion rate from its radiation dominated value -ffrad is to change the depletion factor 
T> by the factor ffrad/-f^b- Increasing the expansion rate decreases the depletion because the sphaleron rate 
decouples at a higher temperature. 

Is this change significant? For a given theory (with all parameters determined) the depletion factor is 
(in principle) calculable. There is essentially no depletion for any expansion rate greater than the expansion 
rate -ffsph given by setting 2? = 1 in ([2^). For H < -ffgph, however, a baryon asymmetry produced at the 
first order phase transition is attenuated by a factor e~^=p'^/^. Whether a change in the expansion rate 
from that in a radiation dominated universe to a different value is important therefore depends on what 
the critical expansion rate -ffsph is in the particular model. If a model has -ffsph = 10"-ffrad) the baryon 
asymmetry left behind in the universe with H ~ .ffsph may be compatible with observation, and that 
in a radiation dominated universe too small by a factor e~^^'\ If, on the other hand, -ffgph < -f^rad the 
asymmetry will survive unattenuated in either universe. 

We now turn to determining the effect of treating the expansion rate as a variable on the sphaleron 
bound in its more familiar forms, in which the requirement P < 1 is converted to a bound on parameters 
in a particular model. The bound is usually stated as a lower bound on the sphaleron energy, or as a 
lower bound on the ratio of a VEV to the temperature at the nucleation or critical temperature, and 
then converted into a bound on parameters in the specific model concerned. We will follow through the 
derivation of such bounds in detail, particularly because we wish to note certain points which are often 
overlooked in this context. We then analyse the case of the MSM in detail using the same effective potential 
(0) and (^ used in the previous section. 

Using the sphaleron energy x = B{Atx / g){cj){T) /T) as the variable in (|2^) we obtain the sphaleron bound 
in its new form as 



sph 



All J 
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di<p/T) 



dT 



(x) 



-1 



dxx'^e 



(23) 
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where, to derive the latter expression, we assumed that over the range of temperatures which contribute to 
the integral the derivative term is approximately constant, and {T\^/T{x)Y ~ 1- Let us assess the validity 
of this approximation in more detail in the case of the MSM. At any temperature at which the two phases 
coexist, i.e. between and Tq, 



T— (t 
dT \T 



where TdXx/dT = — 8/3, and 



27(ro/r)^ 



T dXx 



2Ar(0/r) - a Xt dT 
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(24) 



«+[a2-4AT7(l-(To/T)2)]i/2 



2At 



(25) 



For simplicity, in ( p4D we neglected the temperature dependence of a and 7, which would result in numer- 
ically irrelevant corrections. Even though both (j)/T and its derivative in (24) and ( p5| ) are very sensitive 
functions of T, in the temperature interval Tc > T > Tq both are monotonically decreasing, and we can 
write their lower and upper bounds as follows: 
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(26) 
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where ^ ~ 18. The large value of the derivative means that the pre-factor in front of the sphaleron rate 
in ( p3|) is > 10^. This is essentially just the (inverse) fraction of an expansion time in which the sphaleron 
freezes out (leading to the difference from the naive freeze-out estimate -ffgph ~ rsph(7b)). The range of 
temperatures which contributes in the integral is therefore much less than between T^ and Tq ~ O.OOTc, 
and the constant derivative approximation used in evaluating it is indeed very accurate. Further, as Tb 
varies in this range the change in the result associated with the derivative is at most this factor of three. 
In what follows we will keep track of this dependence of the sphaleron bound on Tt,, and quantify it in 
comparison to the other effects on the bound in which we are interested here. 

The sphaleron bound as given in ( p3| ) can be converted, for a given expansion rate H^, into a lower 
bound on the ratio 4'h/Tiy (where (p]^ = (j)(Ti,)). A numerically convenient and instructive way to write the 
lower bound on this quantity is in the implicit form 



Tb B \4tt 



6Nj.an ( 2TTg<f>[T) MrAArot Vrot^sph 
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dTT 
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ln-^ + 71n^ 



(28) 



In figure 2(a) we show the solutions to this equation obtained from an iterative evaluation of (28), for 
the MSM. We have taken Vrot = Stt^, and fit A/'trA/'rot ^ 8Q - 51n(m|^/8m^) The one loop result 

for /Cnuci we took from fH]: /Csph = {7.54,5.64,4.57,3.89,3.74} for mn = {0.4, 0.5, 0.6, 0.8, l}mvF , and 
extrapolated or interpolated for other values rriH G [10, 80]GeV. u>- we took from |^], where it was found 
that uj-/g(t)iT) G [0.4,0.55] for mn £ [10, 80]GeV. We neglected the plasma effects on uj^. Finally, we took 
an = 0.4, andS = {1.52,1.61,1.83,2.10} for mf^/m^ G {0,8x10-3,8x10-2,0.8,8} H, and quadratically 
interpolated for the intermediate values. 

Besides varying the expansion rate over the range shown in the figure, we have taken a wide range 
of Higgs masses and different values for Tb. It is instructive to do this because the sphaleron bound is 
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Figure 2(a). Minimum <t)(T)/T vs. tine expansion rate 
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often stated as a bound on this ratio of VEV to temperature as if this were a model-independent and 
temperature independent statement of it. We see from figure 2(a) that this is very far from being true. For 
a range of Higgs masses from lOGeV to 80GeV the bound on (j)/T decreases by about 20%. That most of 
the dependence comes from the factor which varies non-neghgibly with the Higgs mass, can be verified 
easily^. This is also clear from figure 2(b), which shows the sphaleron bound as a lower bound on the 
sphaleron energy i^^pj^/Tb (where E^^^i = Esph{T\y)) as a function of the expansion rate. There is however 
still a significant mass dependence (approximately 8% over the mass range considered) in the bound stated 
this way. The temperature dependence of the bounds is comparatively smaller - as T increases from Tq 
to Tc, the bound on (j)/T decreases by 3 — 4%, and for Tb G [Tkucb^c] by less than 1%. This dependence 
comes from the derivative of the VEV inside the logarithm, which as we saw above can vary by a factor of 
three over the range from to Tq. In analyzing any particular model in detail the parameter dependence 
of the sphaleron bound stated this way in terms of these quantities should clearly be borne in mind and 
carefully examined. 

The dependence we are primarily interested in here is that seen in both figures 2(a) and 2(b) on the 
expansion rate of the Universe. Both (/>b/2b and E^^-^/T^h show an almost exact linear dependence on the 
logarithm of H which is evident from (pSj). For a small fractional change in the bound on 4>h/Th or E^^^/Ti, 
due to a change in the expansion rate from i^i-ad we have the approximate formula 



Tb 



In 



H 



Tb 



(29) 



which, given (47r/a^) ' ~ 20, agrees well with the numbers read off from the figures. Over five orders of 



^5{<l)/T)/{4>/T) ~ [l -7S-^(a„/47r)^/V(c?!>/r)] \ For m„ = lOGeV to SOGeV, 5B/B 

5{(j)/T)/{<j)/T) ~ 1/6, accounting for most of the dependence on mn on figure 2(a). 



1/8, and hence 
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Figure 2(b). Minimum sphaleron energy vs tine expansion rate 
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magnitude in the expansion rate we see a decrease in the bound on </>b/Tb by about 0.4, or approximately 
0.08 per order of magnitude increase in the expansion rate. 

The usual starting point for analysis of most extensions of the standard model departs from the 
sphaleron bound given as a lower bound on the ratio of the sphaleron energy or the appropriate VEV 
to the temperature, and then converts this to a bound on the parameters of the model. We have noted 
that such a procedure should be considered more carefully as there can in fact be significant model de- 
pendence in the bounds on these quantities. We have derived nevertheless how such bounds are changed 
as a function of the expansion rate of the Universe, and the approximate form ( |29| ) is essentially model 
independent. Using this formula one can therefore turn the usual sphaleron bound for any given model 
into a lower bound on the expansion rate as a function of model parameters, provided one has the correct 
form of the bounds on (j)h/Th (or E^^^/Tij) in the radiation dominated case: For each set of parameter 
values one calculates the value of the given ratio, and then solves using (^9|) for the expansion rate which 
reduces (or increases) the radiation dominated value to the calculated critical value. 

However, the most direct way to calculate the sphaleron bound as a lower bound on the expansion rate 
is simply to evaluate the integral (l23|) directly to find Hgph for each value of the parameters in the theory. 
We have done this for the MSM using the same parameter values and effective potential as above, and for 
the temperatures Tb = TcTnuch^o- The result is shown in figure 2(c), where the sphaleron bound is given 
as a plot of the minimum expansion rate required as a function of the Higgs mass niH- The dependence on 
the temperature seen in the figure is greater than in the bound on since it also enters in relating rriH 
to (/i/r, as shown in figure 2(d). Figure 2(c) shows dramatically how badly the usual sphaleron bound is 
violated in the MSM ^. For no physical Higgs mass is the minimum required expansion rate within orders 
of magnitude of that in a radiation dominated universe. The discrepancy of this result with the early 
sphaleron bounds calculated for the MSM |jl^ is explained by the much larger (now physical) top quark 

^Studies of the two loop effective potential and lattice studies show that the one- loop ring improved effective potential we 
are using underestimate the strength of the phase transition, but not enough to significantly alter the conclusions drawn here. 
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Figure 2(c). Minimum expansion rate vs m^ 
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mass vrit = 175GeV used here. For small ran one can see from (|8|) that the one-loop thermal contribution 
from the top quark dominates Ay, and therefore, from (p5|), that ^jT stops increasing and levels off as 
seen in figure 2(d). The increase in the minimum expansion rate seen in figure 2(c) as mu decreases below 
this value comes simply from the dependence on mu of the sphaleron energy through B (which decreases, 
increasing the sphaleron rate). 

In many extensions of the standard model it has been shown that, in contrast to the MSM, there are 
physically allowed regions of the parameter space where the usual sphaleron bound is satisfied. The way of 
stating the sphaleron bound we have illustrated for the MSM can be easily generalized to any such model. 
Besides being, as we argue in this paper, a more correct way to state the sphaleron bound (given that the 
expansion rate really is an unconstrained parameter), our discussion also shows that it is an instructive 
way to state it, because it quantifies how well or badly the bound is satisfied or violated. If we state the 
sphaleron bound in this way, it is easy to determine the effect on the calculated bounds in any change to 
input parameters (e.g. to any of the pre- factors in the sphaleron rate). 

Having discussed how the sphaleron bound should be restated as a lower bound on the expansion rate, 
let us ask finally what electroweak baryogenesis at a first order phase transition can potentially tell us 
about the expansion rate at that epoch. A priori we do not know what it is and can use baryogenesis 
as a probe. If the correct electroweak theory turns out to be one in which there is a first order phase 
transition which successfully produces exactly the right amount of baryons during the phase transition, we 
would have compelling evidence that the expansion rate is greater than the corresponding critical value 
Hsph- But it can tell us no more. If the model satisfies the 'old' sphaleron bound with the assumption 
of radiation domination, but has -ffsph < -f^rad (as it typically will), the success of the model provides no 
evidence that the Universe expands at -ffrad- It could even potentially expand orders of magnitude slower 
than -ffrad- We will now see that in contrast electroweak baryogenesis in a homogeneous universe provides 
a much more sensitive probe of the expansion rate at that scale. 
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Figure 2(d). (j)(T)/T vs itih in MSM 
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4 Baryogenesis in a Homogeneous Universe 

Analysis using the effective potential constructed in perturbation theory indicates a first order phase 
transition but is only of validity for Higgs masses up to about 60GeV. Recent non-perturbative results 
indicate that for heavier Higgs masses the line of first order phase transitions ends in a second order phase 
transition at about 80GeV in the Minimal Standard Model. For larger masses the transition is an analytic 
cross-over, i.e. there is actually no phase transition since all physical quantities vary continuously (and 
differentiably) as a function of temperature. This sort of behavior is typical of a system in which there 
is no order parameter which can define the symmetry state of the system - the gauge symmetry is never 
strictly speaking broken or unbroken. 

The only departure from equilibrium in this case is that caused directly by the expansion of the 
Universe. All physical quantities vary on a timescale characterized by the cooling rate ~ H. Unlike the 
case of bubble nucleation there is no separation between the mechanism by which the baryons are created 
and the part of the calculation involving the expansion rate directly, a separation which allowed us to take 
the created asymmetry simply as an input without specifying how it was created. Here we must make 
use of a particular model in order to answer the question of how the baryon asymmetry depends on the 
expansion rate. 

Most work on mechanisms for electroweak baryogenesis has considered extensions of the Standard 
Model with an additional source of CP violation beyond the KM matrix. On bubble walls formed at a first 
order phase transition the CP violation produces in various ways a term biasing the anomalous sphaleron 
processes, causing the creation of baryons on or around the wall. These source terms, which are present 
when there is space or time dependence of the condensate fields, can equally be used to bias the anomalous 
processes and produce baryons when the phase transition is not first order. In the case of a second-order 
or cross-over transition we expect the evolution to be homogeneous with time dependence only of the 
condensate fields, and we will model the problem this way. In fact the validity of the analysis is broader 
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than just the regime where the phase transition is not first order. It also describes well the period after 
the completion of a first order phase transition. In particular, as we will discuss below, it describes the 
case where the phase transition is too weakly first order to satisfy the sphaleron bound. 

We will now consider separately two types of source terms for baryogenesis discussed in the literature. 



4.1 Potentials for Baryon Number 

The first apparently viable mechanisms for electroweak baryogenesis, discussed in p5| and |^6|, considered 
potentials for baryon number. The models differ in their particular realisations of this potential. In various 
theories - two doublet extensions of the minimal standard model |25| and supersymmetric theories with or 



without an additional singlet |^6|] - there are CP odd terms in the effective action for the gauge-Higgs sector, 
of the form {g^ / \Q'k'^)xF F , where F and F are the SU{2) field strength tensor and its dual, x is some field 
or combination of fields which acquire VEVs at the phase transition, times a numerical factor (typically a 
suppression). When these terms are integrated by parts and the anomaly equation {g'^ /167r'^)FF = d^jj^ is 
used, in the homogeneous case (with time dependence only) they produce terms calculationally equivalent 
to a chemical potential for baryon number xB- Specifically in two doublet models there are terms with 



[25] 



XB = -^(^)'4z(fIPo$l-(Po<5l)tcI>i 



ttTJ 

- 7C3f^)'^^^, C3- 1.202, (30) 
\tt1 J vf + vf 

where 6 is the relative phase between the two doublets, with VEVs of magnitude vi and V2 (where the 
former couples to the top quarks). In theories with CP violation characterized by some scale M ||2^ the 
equivalent quantity xb is 

^ rdo\ct>\\ ^dos, (31) 



3M2 ' 3M 

where the first case is a theory with doublets only, the second one with a singlet s. 

Up to higher derivative corrections to the VEVs the system in this background tries to thermalize to 
the equilibrium in the presence of this extra term, in which the baryon number is given by the expression 
in (p^, with fj,B = and H including the additional term due to the background. 

It is perhaps instructive to note exactly why an equilibrium calculation of this sort can give a baryon 
asymmetry. The Sakharov requirement of a departure from equilibrium is sometimes shown using this type 
of expression by acting with @ =CPT on baryon number as follows 

{B) = Tr \@Be^ee-f^^Q^] = -{B) (32) 

and therefore {B) = 0. The same will hold true if we allow non-zero chemical potentials for charges which 
are CP even and, of course, it will not hold if we impose a chemical potential for CP odd charges like B or 
B — L. In the case we are considering the reason it does not hold is that the additional effective term in 
H is not CPT invariant, as the time varying condensate field violates CPT spontaneously. The underlying 
Hamiltonian is of course CPT invariant, but in the expanding Universe this symmetry is spontaneously 
violated. 

The constraints which are to be imposed are those setting all the charges which are conserved over 
the relevant timescale to zero. Compared to the case of the same source term used to generate a biasing 
of baryon number on the bubble walls during a first order phase transition, there are thus two important 
differences: 
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(i) The 'relevant timescale' on the bubble wall (of thickness L moving with velocity v^) is the wall 
passage time Ljvw, typically ~ 10^/T. Here it is that characterizing the time rate of change of the field; in 
the homogeneous case this will be ~ . Thus the numerous processes (e.5. chirality flipping processes 
of the lighter quarks and leptons), which are effectively inoperative on the bubble wall, are equilibrated 
in the present case. The set of relevant conserved charges is therefore much smaller (and the calculation 
therefore simpler) . For a radiation dominated universe the only conserved charges at the electroweak scale 
are the exactly conserved charges - hypercharge y, electric charge Q, and ^B — Li {Li the lepton number in 
generation i). The charge which is violated slowest is right-handed electron number cr, since it is coupled 
to other species only through its small Yukawa coupling y^, by processes with rate ~ y'^T ~ 10~^^T (times 



a number of order 0.1 — 0.01 - see |27], [28]). Thus for H > lO'^^T we will also need to add br as a 
conserved charge, for H > lO^^T both eji and /i^, etc. 

(ii) On a bubble wall the constraints forcing the conserved charges to zero are appropriate only when 
negligible charge can be transported onto the wall over the relevant time scale. (The charges are conserved 
only globally, not locally on the bubble wall unless this is true |^9|.) This places a condition on the thickness 
of the wall L > D/v^ for the applicability of this simple form of the calculation. This condition follows 
from the requirement that the wall passage time be greater that the diffusion time (L/vw > D/v'^), in 
order for transport to be inefficient. In the present case the Universe is (assumed) homogeneous and the 
global constraints forcing the charges equal to zero are always appropriate. 

Because the electroweak phase transition is not a symmetry breaking phase transition, we cannot 
define an exact criterion for whether the 'broken' or 'unbroken' basis of particle states presents the correct 
description in the equilibrium calculation. A correct calculation would assume neither basis. A simple 
example of such a calculation has been discussed recently in |3C] and a cross-over from one limit to the 
other explicitly shown to occur at mw ^ m£, (the vacuum mass and Debye mass of the gauge bosons 
respectively). Converting this to a constraint on the ratio of the VEV to temperature, it turns out that 
the symmetric phase calculation is a better approximation when the sphaleron freezes out. Thus we will 
calculate here in this approximation, using the 'unbroken' phase classification of the states. In section 4.2 
we will also see that either basis of states gives almost identical results. 

We take the case when e/j = as this will turn out to be appropriate in the regime of expansion 
rates of interest for the generation of the observed BAU^. However, the numerical difference induced by 
this additional constraint will turn out to be insignificant. Expressing the charges in terms of particles 
densities n^, and using the linear approxinicition tT/q, — (X'^/12)/cq, ^Q, where is given in terms of the 
chemical potentials ha for the charges Qa by = <la^^A (where is the Qa charge of the species a) and 
ka is a statistical factor which is equal to 1(2) for fermions (bosons) in the massless approximation]], we 
find 
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(10 + n)iXY + 2/iB + -^jfJ.j 



(33) 



^We will neglect the potentially interesting effect discussed in Incorporating it could lead in certain cases to minor 

changes to the final baryon asymmetry. 

Including the lowest order mass correction to this simple formula results in; Ua = (rV6)/i«(I - (3/2)(m„/7rT)2) for 
fermions, and Ua = (r^/3)/ic«(l — (3/2)(ma/7rT)) for bosons. In this paper we ignore these mass corrections since they appear 
as a subleading correction to the result presented in the text. 
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B = — 
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2/iy + 4/iB + -Sj/Uj 



(34) 



where is the chemical potential for B — Li, and n is the number of Higgs doublets. 
Setting the first three charges equal to zero0 we find 

R ^2 1 36(29 + 6n) 

6 399 + 82n 

Note that Cn = Bo/{T'^fiB) is almost insensitive to n, the number of Higgs doublets, varying only between 
0.436 and 0.439, as n changes from to oo 0. Taking the source term to be as given in (30) or (|3l|), we 



set IJ-B = Xb in ( p3q ) to obtain the baryon asymmetry in the "equilibrium" to which the baryon number 
violating processes will try to bring the plasma in the slowly varying background. To calculate the rate at 
which these slow processes bring the system to this state, we impose a chemical potential /i^ on baryon 
number and include the source term. The baryon number B in this state is then given by ( |35| ) with the 
replacement fJ-B ^ fJ-B + Xb, and therefore 

fJ'B = c;,^ - CnXB^ ■ (36) 



2^2 

Using ([TtI) this gives the rate at which B approaches its "equilibrium" as 



2 Tsph 1 



B = -anTsph{B - CuXbT ) , Tsph = QNp^^, a„ = — , (37) 



where an = 0.382 0.380 as n = ^ oo and Tsph is given in (|19D . 

Before calculating the final baryon asymmetry we consider another treatment of induced source terms 
for baryogenesis in a time dependent background. 

4.2 Potential for Hypercharge 

A different treatment of the biasing of baryon number was given in ||3^|^. In the broken phase of a two Higgs 
doublet model the relative phase 9 of the neutral components of the Higgs fields enters in the fermionic 
mass terms. A hypercharge rotation of the fields to remove this phase from the mass term produces at tree 
level a real mass term and an additional term in the Lagrangian which, in the homogeneous case can be 
written simply as Xy^, where Y is the hypercharge operator |^, and 

2i,2 

XY = -^-^0- (38) 

®In the 'unbroken' phase one can choose to constrain any two hnear combinations of hypercharge Y and isospin T3. The 
choice of Y and T3 is simple, because T3 is then proportional to its own chemical potential and T3 = is trivial. 

^If we had not assumed the right-handed electron to be in equilibrium, the change would be small. In this case Cn = 
0.455 0.462, as n = ^ cxi. 

'^^The account given here is not precisely that of the original version of the idea given in which treated a potential 
for fermionic hypercharge. It was pointed out in that the rotation should also be performed on the Higgs fields, giving 
a potential for total hypercharge which in the unbroken limit (VEVs —> 0) is pure gauge, and therefore can have no physical 
effect. The leading baryon production is in this case mass-squared suppressed. For a discussion of this point, see also section 

We follow the convention used in [p2|. There is nothing special about the choice of hypercharge. The essential element 
is that it is an anomaly free charge which is spontaneously broken by the mass term. A rotation proportional to isospin, for 
example, or any charge which is a linear combination of hypercharge and a charge exactly conserved in the broken phase is 
equally good. It is not difficult to check that the extra induced 'source' term always drops out in the calculations given below. 
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In the unbroken phase this is just a gauge term, but in the broken phase it can have physical significance 
because hypercharge is not conserved, being violated by VEV suppressed terms. 

Again, as discussed in section 4.1, we can calculate either in the broken phase or unbroken phase, but the 
latter is probably more appropriate for the temperature range of relevance. In this case of course we must 
include the information about hypercharge violating processes to get a non-zero answer, so calculating in the 
unbroken phase means taking the basis of chiral states of the unbroken phase and treating the mass terms as 
interaction vertices which can violate hypercharge by flipping chirality. The correct constraint calculation 
is therefore one in which we take the same global conserved charges as in the previous calculation, but 
instead of the constraint on hypercharge, we must impose the constraint on the (conserved) electric charge 
(5, and we get 
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(39) 



The "chemical potential" /iy for hypercharge here is the effective one which arises from the source term 
for hypercharge, i.e. jiy = —Xy- Setting the conserved charges to zero in (^9[) we can solve for B to find 
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Using ( |18|) we obtain the equation describing the relaxation of baryon number to the "equilibrium" in 
presence of the hypercharge source term: 
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(41) 



Note that c'^ = 0.033 0.024, and a'^ = 0.3805 0.3796 for n = ^ oo. The rate of relaxation is 
essentially independent of the number of doublets, and it is almost identical to the rate in the presence of 
a potential for baryon number in section 4.1. 

It is noteworthy that the coefficient is significantly smaller than the corresponding coefficient in 
the case of a potential for baryon number (c„/c^ ~ 13 — 18 as n = — > oo). It follows therefore that, 
even though the baryon number source (|30| ) is suppressed by a factor of mass over temperature squared 
relative to that in (^8|), the former may give the dominant source term for baryogenesis. The reason for 
this is a suppression due to strong sphaleron processes in the case of a hypercharge source term [35|. In 
the massless quark approximation these force the densities of right and left-handed baryons equal, i.e. 
Bl = Bpi. On the other hand, it is easy to show that, with source terms for a charge such as hypercharge 
which is conserved in baryon number violating processes, B oc {3Bl + Ll). With B — L = this implies 
B oc (|i? — L/j) (where L/j is the density of all right-handed leptons). Therefore, setting Lji = would 
lead to an equilibrium with B = 0, i.e. a vanishing source term for baryon number. The non-zero result 
we have obtained is therefore proportional to the charge in the right-handed leptons, which gives a small 
statistical factor related to the fraction of the total number of degrees of freedom they represent. 
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When the VEVs approach zero, the result in ( pT[ ) does indeed vanish, but not explicitly. In this limit 
the rate of the 'hypercharge violating' processes goes to zero, so the "equilibrium" calculation is no longer 
appropriate as it is only valid on a time scale longer than one which diverges as one sends the VEVs to 
zero. It is not difficult to check that one does indeed get zero for the baryon number in the presence of 
this source term when we impose the extra constraint appropriate in this limit |^. 



4.3 The Baryon Asymmetry 



The equations p7| ) and 
together. Integrating (| 

B{t)- 



|4l] ) have the same form for both source terms, and so we can analyse them 
gives the baryon asymmetry as a function of time as 



dt"anTsph{t'' 
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(42) 



where tj is an initial time chosen before the phase transition or cross-over takes place such that the source 
term may be taken to be zero. This expression is simply the source term integrated against the appropriate 
Green's function. The freeze-out time (temperature) tf(Tf) is that at which the integral in the exponent 
is equal to one, i.e. 

/* dt"a„rsph(t") = 1 . (43) 
Jtf 



The approximation in (42) follows since we would expect that the time-scale characterizing the variation in 
X should be of the same order as that characterizing the change in (j)/T. However, as discussed in section 
3.2, there is an exponential dependence in the sphaleron rate on (p/T with a large pre-factor (~ 1/aw). 
This means that the derivative inside the integral in (^2|) can be approximated by a delta- function at tj 
down to a time scale much shorter than that over which x varies, and the result follows. The sphaleron 
rate T^ph enters only in determining the freeze-out value for the source [c„x^^]frccze- Optimally, the 
sphaleron processes switch off when the source is at its maximum, leading to an estimate of the maximum 
production of baryons at a second order phase transition or cross-over: -^^rnax ~ [^nXT ]max* To a very 
good approximation the final baryon to entropy ratio is 
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(44) 



using the fact that dx/dt = —HTdx/dT and that the entropy density s = (27r^/45)g^,r'^ {g^, the number 
of relativistic degrees of freedom). The subscript denotes that all these quantities are to be evaluated at 
the sphaleron 'freeze-out'. 

In the case of baryon production in a homogeneous Universe with source terms of this type the final 
baryon asymmetry is therefore proportional to the expansion rate at freeze-out. This contrasts completely 
with the case of baryogenesis at a first order transition, for which the baryon asymmetry can be effectively 
the same for an expansion rate differing by many orders of magnitude. 

We can invert (|4j) to get the range of expansion rates consistent with the baryon to entropy required 
by nucleosynthesis^^ 

~(2-12) — — J— r (45) 



freeze 



\{Tdx/dT) 



freeze I 



Choosing to impose the conservation on T3 with chemical potential ^Ts, we get T3 = -^(10 + n){fiT3 + ^q) and Q also 
picks up the extra term (10 + 71)^X3- Imposing T3 — leaves only the linear combination 9 + fig in the other equations. The 
solution is the trivial unperturbed equilibrium. 

""^^This range corresponds to the conservative bounds from direct observations of element abundances given in [^. Tighter 
bounds, corresponding roughly to the range (3 — 9) x 10~^^ in ([45|), are given in and ]3q].. 
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How big is {Tdx/dT)fT.ceze in any given theory? A full treatment of the phase transition in any of the 
models mentioned would be required to actually calculate this, a task however considerably beyond the 
methods used to date in the study of the phase transition^. A naive guess would be 

~ ^6 ~ e (46) 

taking the field x to trace the VEV (or combination of VEVs), which is itself then assumed to evolve roughly 
in proportion to the temperature {i.e. on a timescale given by the expansion rate). The parameter e is one 
characterizing CP violation, which we might expect to be constrained by CP violation phenomenology of 
the relevant model. 



A full calculation of any given model at finite temperature would be required to turn the bound (45) 
into a precise one on the expansion rate alone. However, short of such a calculation, we can do better than 
the very naive estimate given by (^6|). 

(i) In section 3.2 we examined the Minimal Standard Model and saw that, near the critical temperature, 
the VEV is a very sensitive function of the temperature, with Td{(j)/T) / dT ~ d(f)/dT ^ (100 — 30) in the 
range of temperatures to Tq, and about 60 at the nucleation temperature. Typically the sphaleron will 
freeze-out in this range of temperatures, as the sphaleron rate changes by many orders of magnitude. The 
same sort of behavior can be seen to continue at larger Higgs masses in the non-perturbative treatment of 
the phase transition, in the case when the phase transition is a "sharp - but regular - cross-over" |^^. This 
means that the range of temperatures over which physical measurables like the susceptibility vary is a small 
fraction of the temperature at which the change occurs. (It is, of course, this "sharpness" which allows 
one still to talk about a phase transition when, strictly speaking, there is none.) From the data in ||2^ we 
see that there is a range of temperature of a few GeV which compares with a "transition temperature" 
anywhere between 60 and 200GeV. Thus the standard model estimate of d(j)/dT ~ 60 seems reasonable, 
much larger than our naive estimate in (|46|). 

(ii) We can also learn something about the constraints on e by looking at the effective potential for a 
particular case. Consider a two Higgs doublet model. One interesting regime is that in which the evolution 
of the CP violating angle is determined dominantly by terms which break CP spontaneously. In a CP 
invariant Higgs potential [|39[| only the terms 



\^{{<p\(t>2f + h.c.) {\Q<l,% + \j(t)\<l}2){<t^% + h.c.) (47) 

are functions of the relative angle of the two VEVs. Taking the real parts of the VEVs to be determined 
by the rest of the potential (i.e. working in the approximation that the terms ( |47| ) are small) this gives a 
quadratic potential for the cosine of the relative phase, which (taking A5 — A7 positive), is minimized at 

cos 6 = -mm 1, h — ; (48) 

V ^\5V2 4A5fl/ 

There are two (CP conjugate) solutions which will be split by additional explicit CP violation. How the 
angle changes as the VEVs do depends on the values of the ratios of the couplings Ag/As and A7/A5. 
A necessary condition for dO/dT 7^ at the phase transition is cos 9 7^ —1, which is the case when 
A7/2A5 < vi/v2 < 2A5/A6. For couplings such that the first term in (48) dominates, and vi changing faster 
than V2 as a function of temperature, we have 

rj,de ^ 1 T dvi !!l>/^^V d\nvi ^ d\nv2 ^^^^ 

dT tan 9 vi dT ' V2 \ Ag / ' d In T d In T 



^^Perturbative methods apply when the phase transition is fairly strongly first order. The methods which have been employed 
to describe the opposite regime do not include the evolution of the CP odd fields relevant here. 



21 



Typically we have Td9/dT ~ dvi/dT, but there are also parts of parameter space (near cos 6 = — 1) 
where the phase can change much faster than this. The only role of the explicit CP violation here is to 
split the two degenerate minima so that the same sign is chosen everywhere. This illustrates that the 
constraint on the parameter which we called e from CP violating phenomena at zero temperature may 
be extremely weak. With a moderate fine tuning it can be considerably larger than one, and not related 
directly to any small parameter associated with the smallness of CP violation. In fact in theories such as 
the Minimal Supersymmetric Standard Model (MSSM) it is naturally the case that the terms which break 
CP spontaneously (which are induced in the plasma through thermal corrections) are dominant over the 
terms which break CP explicitly (which are suppressed by a loop factor) |4^. We conclude, on the basis 
of a simple analysis of the two Higgs doublet model, that the naive estimate (^) for with e ~ 1 is 
too small by about two orders of magnitude. A result of this magnitude is obtained for a large portion of 
the parameter space, without any tuning. With a moderate fine tuning, the effective CP violation can be 
further enhanced. To make a more precise statement would require a detailed analysis of the Higgs sector 
of the particular model. 

5 Non- Standard Cosmologies 

Having established quantitively the dependence of the baryon asymmetry on the expansion rate in two 
possible scenarios for baryon production at the electroweak scale, we now turn to the discussion of physical 
mechanisms which could lead to such a different expansion rate at the electroweak scale. 

As mentioned in the introduction this kind of question has previously been treated in the context of 
calculations of relic densities of weakly interacting particles in [^] and The relic density of a weakly 
interacting species is determined by the temperature at which the species decouples from the ordinary 
(visible) matter, which depends, just as in the case of the sphaleron decoupling discussed above, on a 
comparison between the appropriate interaction rate and the expansion rate of the Universe. In typical 
models this decoupling occurs before nucleosynthesis, and therefore one is led to consider, just as we are 
doing here, possible alternatives to radiation domination at that epoch. The alternative which is considered 
exclusively in Q and in most detail in Q is: 

(i) An Anisotropic Universe: A universe which is homogeneous but not isotropic is described by a metric 
with three scale factors, one for each spatial dimension. With an adiabatic approximation (i.e. expansion 
slow enough to allow thermalization) it is the effective volume expansion rate H associated with an average 
scale factor a which determines how the temperature changes in the same way as in the isotropic FRW 
spacetime. The equation of motion for H is just that of the FRW space, but with an additional term which 
is equivalent to a component of the energy density scaling as 

A further possibility considered in [Q] is 

(ii) Non-standard theories of gravity. The case studied in Q is a Brans-Dicke theory, which again turns 
out to effectively produce an extra component in the energy density scaling as l/o^. There are also of 
course many other variants on standard Einstein gravity which can be considered. 

The simple possibility we will concentrate on is: 

(iii) Einstein gravity with isotropy and homogeneity, but with an extra contribution to the energy 
density which is important prior to nucleosynthesis. As noted by one of us (MJ) in any mode of a scalar 
field dominated by its kinetic energy has the required property, as its energy density can scale as fast as 
1 /a^ . The electroweak phase transition could potentially occur during a phase of the Universe dominated 
by the kinetic energy of a scalar field, termed kination in which can end before nucleosynthesis as the 
kinetic energy density red-shifts away relative to the radiation. Below we will discuss several ways in which 
such a phase can come about within the context of inflationary cosmology, which explains the assumed 
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isotropy and homogeneity. In particular we will discuss models which come naturally out of an alternative 
to the usual theory of reheating discussed by Spokoiny in [Q]. 

A clear motivation for considering such models follows from the calculations in the previous section. If 
we have such a component scaling as 1/a^, the expansion rate is given by 



a J \ a 



(50) 



where is the scale factor when the density in the mode becomes equal to that in radiation and pe is the 
total energy density at that time. The factor /(a) accounts for the effect of decouplings, and, assuming 
adiabatic decouplings, /(a) = [g{a(.)/g{a)]^^^ , where g{a) is the number of relativistic degrees of freedom 
as a function of the scale factor a. Nucleosynthesis constraints place a lower bound on Te, the temperature 
at the time of equality of radiation-kinetic energy density, which can be inferred from the corresponding 
bounds on additional relativistic particle degrees of freedom. This is the case since the predominant effect 
of such extra degrees of freedom is also in the change they cause to the expansion rate at the beginning of 
nucleosynthesis, which determines the crucial ratio of neutrons to protons when the weak interactions drop 
out of equilibrium at ~ IMeV. We take here the conservative bounds of |36], which allow the equivalent of 
1.5 extra Dirac neutrino degrees of freedom over the three degrees of freedom of the Standard Model, i.e. 
we allow an additional energy density at IMeV which is 3(7/8)/10.75 ~ 0.25 of the standard model one 
(with 10.75 degrees of freedom). From (]50| ) this means 



< 0.5 



(51) 



Using Ta = f{a)T(,ae this gives the upper bound on the expansion rate at the electroweak scale 



-< 1.8x10-11 
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(52) 



The result differs only by ^3/2 if we take the less conservative nucleosynthesis bound of |37| and 

Taking i/max to be the expansion rate corresponding to the upper bound (52), the requirement (45) 
for generation of the observed BAU at the electroweak scale in a homogeneous Universe can be expressed 
as a requirement of the relevant CP violating parameter 
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Absorbing the nucleosynthesis limit, i.e. taking H = -ffmax we have the strict lower bound |^ 
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(54) 



From the analysis in section 4.3 it follows that this bound may indeed be satisfied in extensions of the 
MSM such as those we have discussed, without any fine tuning. If the upper bound ( |52|) on the expansion 
rate is not saturated, the CP violation parameter is required to be larger as given by (|53|). As discussed in 
section 4.3, with some fine-tuning of parameters in the potential, this parameter can indeed be enhanced to 



Note that taking the upper bound on the expansion rate in ( |52[ ) corresponds to absorbing the upper bound on effective 
number of degrees of freedom at nucleosynthesis, which is only consistent with the lowest baryon to entropy ratio (increasing 
the expansion rate increases the fraction of Helium) . 
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considerably greater than the typical value which just satisfies the lower bound (54). An exact statement 
of how large it can be would require a detailed examination of the model in question. 

The important result is that in a cosmology with an additional component scaling as 1/a^ which dom- 
inates prior to nucleosynthesis, the creation of the baryon asymmetry is possible (i.e. consistent with all 
observations) at the electroweak scale in a homogeneous expanding Universe. The fact that generation of 
the BAU in this case has generally been dismissed as impossible provides clear motivation for the consid- 
eration of such cosmologies in greater detail. Certainly also as experiment pushes the bounds on scalar 
particles upwards, the usual sphaleron bound for generation of the BAU is becoming increasingly severe 
and alternative mechanisms for the production of the BAU within the context of electroweak cosmology 
become more relevant. 



5.1 Kination 

Consider the dynamics of a real scalar field (p with potential V{(j)). Varying the action 



d x^—g 



(55) 



and taking the FRW metric ds^ = dt^ — a{t)'^dx'^ with scale factor a{t), gives the equation of motion for 
the homogeneous modes which, after multiplication by 0, can be written as follows 
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Defining ({t) = V {(j)) / p{4>) , where p{t) = 4? + V{4>), we find 

/■* t"- da 

p{t) = p{to) exp - / 6[1 - at)]H{t)dt = p{to) exp - / 6[1 - C(a)]— . 
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When the kinetic energy dominates, C ~^ and 
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If a potential possesses a flat direction, for example, the energy in the associated coherent goldstone mode 
scales in this way. In this case the scaling can be seen to follow directly from the conservation of the 
Noether current associated with the symmetry. Consider for example a complex scalar <!• with a potential 
invariant under the global symmetry $ — > e*^$ 
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The mode $ 
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with 9 = const is a solution of the equations of motion for which the conserved 



Noether charge is 

f = pe = a^i^^ (9° $ = 2a^v'^e . (60) 

Thus 9(xl/a^ and p = ocl/a^. 

Such kinetic energy dominated modes represent the opposite limit to inflation [^] which is driven by 
potential energy so that C ~^ 1 ^"^^ p{t) ^ Pi'to)- Indeed for any homogeneous mode (assuming only that 
V{(j)) is positive) we have that 
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Instead of superluminal expansion in inflation, a kinetic energy dominated mode of a scalar potential drives 
a subluminal expansion very similar to that of radiation (a oc t^/"^) or matter (a oc t^/^). Writing the stress 
energy tensor in terms of the pressure p and the energy density p in the standard way, the equation of state 
\s p = p for the kinetic mode, in contrast to p = (l/3)/9 (radiation), p = (matter), and p = —p (inflation). 

We now consider various ways in which a phase of kination could come about. Inflation is the standard 
paradigm which explains isotropy and homogeneity of the Universe as it appears today. A scalar field 
drives a period of inflation and subsequently decays, filling the Universe with radiation and matter. We 
will assume that a period of inflation produces the isotropic and homogeneous Universe, but ask how it 
might come about that after inflation a reheated universe would be dominated by a kinetic scalar mode. 

Two questions can be separated: 

• How can a scalar field potential support a mode that is kinetic energy dominated? This question is 
twofold. Firstly, what shape must the potential have to keep the kinetic energy dominant? Secondly, what 
is required of the field in order that energy does not leak out of the coherent mode? 

• How can kinetic modes come to dominate the energy density, i.e. how can they be excited? 

Let us start with the first question. The most trivial case of potential energy domination is the example 
used above of an exactly flat potential. This case is not of interest here, since the energy in such a mode 
is negligible at the end of inflation, as it also red-shifts away as during inflation. Any kinetic energy 
domination must therefore occur through the roll of a fleld in a potential after inflation. The dynamics of 
a homogeneous real scalar field cf) with potential F(0) in an expanding FRW universe are described by the 
equations 

^ + 3i/<^ + y (<^) = J_^(a3</.)+y'((/>) = (62) 
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Pr+^Hpr = (64) 

where pr is the energy density in radiation, to which we assume the scalar field to be coupled only through 
gravity. This is a roll damped by the expansion of the Universe and the first question is therefore: how 
steep must a potential be in order that the roll be more effective in creating kinetic energy than the 
damping is in attenuating it? A hint of the solution is immediately given by considering again the trivial 
case V{(f)) = 0, which gives the solution 

m = 4>o ' = 0o (y) , m =<i^o + ^oto In ^ (65) 

when pr = 0. If, with this solution for (j), the potential is such that the terms which depend on it decrease 
faster than the other terms in the equations of motion, the kinetic energy domination will continue once 
established. Given that the time dependence is logarithmic, it is clear that an exponential potential is what 
is required. That exponential potentials define what is steep enough for kination can be seen by taking 
the explicit attractor solutions given in El for the potential V{4>) = Voe'^'t' 



2 2 

(j){t) = -\nt, a(xt^, C = l , A2<6, (66) 

A 6 

and the origin of (j) is redefined so that Vq = 2X~^ (6A~'^ — l). (For simplicity we took Mp = 1.) The 
context within which [^] discussed this potential was "power-law inflation", for which the superluminal 
expansion occurs when A < \/2. From (57) it follows that p oc a^'^ so in the limit A \/6 we recover the 



25 



scaling of an exactly flat potential. When A > \/6 the pre-factor cannot be written in this way. There is 
no single attractor solution, and the ratio ^ rather than being fixed approaches zero asymptotically. 

In this analysis we have assumed a simple roll down a potential. Another possibility is that a field 
oscillates about a minimum. It is easy to see from (^) the well-known result that an oscillation in a 
quadratic potential gives an energy scaling like matter oc l/o^ since C can be replaced by its average 
((") = 1/2 over a time-scale of the expansion time. The analysis for the cp'^ potential is given in [43|. The 



result is that an oscillating mode scales as a"^"/^"^^-* and, correspondingly, {0 = 2/(n + 2), so that the 
kinetic energy becomes more dominant as n increases. 

We also require that, if such a mode is excited, the energy remains in it, i.e. that it does not leak out 
by decay of the coherent mode into particle excitations of itself or other fields to which it is coupled. In 
the present context we want to maximize the effect of the mode and therefore need the energy to stay in 
the kinetic mode from before the electroweak scale until shortly before nucleosynthesis. Potentials which 
support coherent modes which are so weakly coupled to other fields (or self-coupled in the case of massless 
fields) that they do not decay before nucleosynthesis are in fact commonplace in particle physics - they 
are the source of problems like the Polonyi problem. In particular, exponential potentials which arise in 
theories involving compactifications are typically extremely weakly coupled to other sectors. Accordingly, 
these sectors are termed "hidden," as they interact with the "visible" matter only (or predominantly) 
gravitationally. We will thus assume that perturbative decay is negligible. On the other hand non- 
perturbative decay mechanisms like parametric resonance which have been much discussed recently in the 



context of the problem of reheating after inflation |5^, 54, 53, 52, must be considered in the case 
of oscillating modes of non-linear potentials. This is a possibility we will consider explicitly below, and 
the requirement that such a mode survive until nucleosynthesis will place constraints on the potentials we 
consider. 

The second question above concerned how such a mode would come to dominate over radiation. In the 
analysis just given of a field rolling down a potential we set pr = 0, and the solutions are therefore valid 
only if » Pr (where prj, is the total energy in the scalar field). What happens if pr ^ /O0? Can we end 
up rolling into the kinetic energy dominated mode with ^> pr if we are in an exponential with such a 
mode? In a radiation dominated universe a cx t^/^ {H = l/2t) so that the damping is stronger than in 
the kinetic energy dominated universe where a oc 

ti/3 (^H = l/3t). If the scalar field is in a kinetic energy 
dominated mode, we find, assuming radiation domination, that 

3 

(^(t) = M^f = (y) ' , m =^o + ikto (l - (y)^) . (67) 

The result is that the exponential potential energy will always "catch up" with the kinetic energy and 
the field will be driven into a mode which scales much slower than radiation, until p^ ~ pr- These quite 
different behaviours in the two limits p^ » pr {P(j) oc 1/a^) and p^ Pr {p ^ const) which tend to drive 
the system from one regime to the other suggests that there may be an attractor solution with ~ 
That such a solution exists and is an attractor has been noted in 113] and |45]. It can in fact be generalized 



to the case that the non-scalar component scales as pm oc 1/a™ {e.g. non-relativistic matter with m = 3), 
and all components scale as it does with the ratios of their contributions given by 

6 P4> + Pm y 

In this case one need not assume that the pre-factor in the exponential can be written in the special form 
required for the solution (^) and A can take on any value A > \/m, which is just the requirement that the 
attractor mode in (B^) with pm = scale faster than l/a^. 
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The existence of this attractor means that if we start in a radiation dominated universe (or, more 
generally, in a universe dominated by energy scaling as 1/a™) we will always end up in this solution ( |68|) 
rather than in the kinetic energy dominated mode of the exponential. In order to realize kination in this 
potential we must therefore satisfy the condition at the end of inflation, or some time after it, that the 
kinetic energy dominate over the radiation. The dynamics of the simple exponential alone will not produce 
kination if we have a radiation dominated universe after inflation. We will examine two possibilities: (i) A 
non-standard theory of re-heating in which the radiation in the Universe is that created by the expansion of 
the background during inflation, and radiation is naturally subdominant after inflation, and (ii) standard 
reheating with a slightly different exponential potential which can first cool the radiation with a short 
inflationary phase and then roll into a kinetic mode. 

We will concentrate on this first kind of model, because in it there must by construction be a phase of 
kination, and in our view it offers a very attractive (and unjustifiably neglected) alternative to the standard 
re-heating scenario. What sort of model would lead to this phase ending as late as nucleosynthesis will be 
the question which interests us in the specific context of electroweak cosmology. 

In the oscillatory potential things are slightly different. The scaling was predicated on the assumption 
that the field oscillated on a time-scale short compared to the expansion time, but not on any assumption 
about the time dependence of the expansion rate (i.e. about which component dominates the energy 
density). Thus if the Universe is radiation dominated when we enter the oscillatory mode of a potential 
with n > 4, it will always be radiation dominated since the energy in the scalar field red-shifts away 
faster. To realize kination in this potential we therefore require the radiation to be sub-dominant when 
the oscillatory phase begins. Just as for the exponential potential we will discuss in this case how this 
condition can be realized in two ways: (i) in the same alternative standard theory of reheating after a 
period of inflation driven by the power-law potential itself, and (ii) with ordinary reheating by another 
field followed by a brief subsequent period of inflation as the field with the power-law potential rolls before 
it begins oscillating. Again most of our attention will focus on the first case, in which a single field is both 
inflaton and 'kinaton'. 



5.2 Reheating by Kination 

Reheating after inflation is required in order to match the 'cold empty' Universe left behind by inflation 
onto the radiation dominated one which must in the standard cosmology be established by nucleosynthesis 
at the latest, (and usually, it is assumed, at some temperature high enough to support some theory of 
baryogenesis) . In the standard theory this is achieved by the decay of the inflaton into particles in an 
oscillatory phase after inflation, the zero entropy coherent state producing the enormous entropy of the 
radiation dominated universe. That there exists a simple alternative to this scenario has been pointed out 
by Spokoiny in [|^]. The Universe is not in fact in an exactly cold zero entropy state after inflation - besides 
the energy in the inflaton, there is some energy in the particles created by the accelerated expansion. The 
process which gives rise to the perturbations from homogeneity required for structure formation on large 
scales creates an energy density which is peaked at the scale H, where H is the expansion rate during 
inflation, with energy density 5p^ = €pH^, where the superscript H denotes that this energy density is 
dominated by the scale k H , and ep ~ {i^'^gt^^ /ZQ)/{2t^Y [HIj where g'f^ ^ 10^ is the effective number 
of light (mass m < H) degrees of freedomf^. In a typical inflationary model with Gncrgy d.Gnsity p\ ^ q-^(^ at 
the end of inflation 



^Pi^,end ^ ^i, end ^ fp Pi , end 
A, end "^3 M2 ~ 9 



(69) 



^This estimate assumes the same contribution from all particles as from the scalar particles analysed in 
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(70) 



which is very small since the energy scale associated with inflation is typically required to be well below 
Mp, e.g. for chaotic inflationary model in a potential A(/>^, , end ~ ep^^i while the requirement 

that one gets density perturbations of the correct magnitude (on COBE scales dp/p ^ 5 x 10~^) gives 
A ~ 10~^^. In the context of ordinary re-heating this small fraction is irrelevant as it is swamped by the 
radiation created by inflaton decay. The possibility envisaged in ^ is one which is easy to see given the 
observations of the preceeding section on kination: If, instead of decaying, the inflaton rolls into a potential 
in which its energy density scales as l/a"^ with s > 4, the energy density Sp^ will come to dominate at some 
time t-k , end after inflation when the scale factor has evolved to ak , end from a; ^ end at the end of inflation 
with 

ak,end ^ /9_JVf^\^ ^ (l]^' ( Mp 

di , end \ Pi , end / J \ ' '^^'^ J ' 

where ak,end is the scale factor at the end of kination (during which p oc l/o*), the phase which inter- 
polates between inflation and radiation domination. The energy in the inflaton simply red-shifts away 
instead of decaying. As discussed in in order to accommodate nucleosynthesis there are two require- 
ments which must be fulfilled: (i) the radiation must thermalize at a temperature above IMeV and, (ii) 
the transition to radiation dominance must occur sufficiently long before nucleosynthesis to satisfy the 
appropriate constraints at that time on the expansion rate. Taking kcs{a) to be the typical energy of 
the created radiation as a function of scale factor, we have k^g = end^^i, end/^J- Assuming that the 
dominant form of this radiation is in standard model degrees of freedom, the interaction rate for processes 
coupling them is ~ a'^k^g (for kes S> M\y and a ~ 1/30 — 1/50). Comparing this to the expansion rate 
« 2/st -ffi, end(oi, end/o)'*''^ (in kination), we get an estimate for the thermalization temperature 

^reheat ^ 

^reheat ~ Hi end ( o ^^^^ ) a~ (71) 

9*,i, end / 

where ff*, i, end is the number degrees of freedom which are relativistic at k ^ H^^ end? and we have defined 
^Pi^ end ~ ff* , i , end'^T^^'^ end/^O, and taken T OC l/o f^. Assuming this temperature Tj-eheat to be attained 
before the transition to radiation dominance, it follows from ( |70D that T^^end) the temperature at the 
beginning of radiation domination, is given approximately by 



1 / \ i 3-2 

2k, end f£p\''-^ I _30ep \ f Hi^cnd^ 



Mp V 3 / \ ir^g^ , i , end / V Mp 



(72) 



Requiring this to be above the nucleosynthesis temperature IMeV places a lower bound on -ff; end- 

For 

s = 6 we find that -ffi^end > 10'' GeV, which corresponds to reheat > lO^GeV, consistent with the 
assumption that Tj-eheat > ^k.end- For s = 5 both -ffi^end and Treheat are greater by a factor of ~ 10^. In 
both cases a late transition to radiation dominance implies that the energy scale at the end of inflation 
and thermalization scale are well below the GUT scale. 

In standard inflationary models the usual constraint on Hi ^ end or the energy density at the end of 
inflation comes from the requirement that the amplitude of perturbations be that required for structure 
formation. In the models which we discuss below we will consider how this non-trivial constraint is 
satisfied in this model of reheating (a question not considered in |||]), and in particular how it fits with 
the particular type of realization of this model we are interested in, where the transition to radiation 

'^^A far from equilibrium system may in fact need many rescatterings {i.e. Nucatt 2> 1) to fully thermalize. Modifying the 
estimate in ^r\\) to incorporate this gives Trohoat smaller by a factor N~'^ll;''~^K 

^^Here and below we neglect the effect of possible decouplings between Trchoat and Tk, end, i.e. we assume the number of 
relativistic degrees of freedom to be fixed. 
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domination does actually occur close to nucleosynthesis with the potentially important consequences for 
electroweak baryogenesis discussed in the first part of this paper. 



5.3 Infiation-Kination in an Exponential Potential 

As discussed in section 5.1, a simple exponential which gives rise to the kinetic energy dominated mode 
required for kination does not itself accommodate an inflationary solution. We need to have a potential 
which is flatter in some region (for inflation) and sufficiently steep (for kination) in the part of the potential 
the field rolls into after inflation. An example is an exponential ~ e~^'^ where A varies as a function of (f). 
As a simple case of this, which we can treat analytically, we consider]^ 

V{(j)) = Voe-^'*' , where A < ^2 for < end 

A = A'>2 for 4>>(t>i,end (73) 

where we set Mp = 1. As discussed above, one solution to the equations of motion for this potential is a 
power-law inflationary attractor (66) with (f> = jln ' ^"-^ ^ 

t^A'. We assume the field to evolve 



in this attractor in inflation from (p <^ (pi^ end • When the field reaches (pi ^ end inflation ends and after a 
transient period it will roll, for A' < \/6, into the new attractor in the steeper potential. If A' > -v/6, there 
is no single attractor, but the field will run after a few expansion times into a solution in which the kinetic 
energy is very dominant. In either case we will neglect the details of the few expansion times in which this 
transition occurs. 

We calculate first the cosmological perturbations generated in the usual way by the amplification 
of quantum fluctuations during inflation. The amplitude of the perturbation in a mode with comoving 
momentum k when it re-enters the horizon after inflation at time 2X, is given by the usual formula 

'Ak)^esf^, (74) 

p 4>ix 

where IX denotes the time when the perturbation k exits the horizon in inflation, and = 3/57r ~ 0.2 
is a constant. The formula is valid provided the slow-roll approximation holds at this time, which in the 
case of the exponential potential corresponds to A <C \/2- The evolution of scales is illustrated on figure 
3, where k = A;physicai« is plotted versus Ha (both on the logarithmic scale). Since the comoving scale is 
fixed, it follows that 

k = kix = k2x <^ {Ha)ix = {Ha)2x , (75) 

and therefore 

{Ha)ix (ffa)k, beg {Ha)r, beg (-ga)m, beg _ ^^g^ 
(-f^a)i,end (-f^a)k,end (-f^a)r , end {Ha)2X 

where the indices i, k, r, and m denote inflation, kination, radiation and matter, respectively, and we 
have assumed that {Ha)\^ end = {Ha\^ beg, [Ha\^ end = {Ha)r ^ beg, {Ha)r ^ end = {Ha)^, beg, and that the 
relevant perturbation enters the horizon in the matter era. In writing ( |76|) we assumed sharp transitions 
z ^ A; — > r — > m. Within our approximation we keep a continuous, but its derivative exhibits a jump 
{H = 2/A^t in inflation matches onto H = 2/A'^t in kination). With the attractor solutions for (^) in 
(|7^) we obtain 

— (^) = e^Ti:;— • (77) 
P A'^ tlx 

^^"We could of course consider any potential which accommodates inflation in some region and is asymptotically a sufficiently 
steep exponential. Motivation for an exponential with varying A is given in [M. 
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Figure 3. Evolution of scales in the inflation-kination model 
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Using ([76D and calculating Ha in each of the eras, we find that 
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(78) 



where, for simplicity, we take 2X to be in the radiation era. The behavior of (^) is the usual one, with an 
overall amplitude set by the expansion rate at the end of inflation and, for sufficiently small A (required for 
consistency of our slow-roll approximation), a fairly flat spectrum of perturbations over the scales relevant 
to structure formation. Indeed, using the standard expression n « 1 + [—3{V'/V)'^ + 2V" /V] Mp |46|, we 
obtain n « 1 — A^. Assuming the reheating scenario of M discussed in the previous section, we can express 



(|77| ) in terms of the radiation temperature at the end of kination, end- Taking the radiation energy 
density at the inflation-kination transition to be Pi,end = ^pH^ end' ^^'^ using (70) and (72) allows us to 
express the temperature and time at the end of inflation as follows 



pA'2-2 
" i , end 
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i , end^p 



^ k , end ' 



Hi. 



end 



^ 9* , i , end 

30Z 



end 



A^ti. 



end 



(79) 



so that 




A' ^- 



-'k,end -^2X 



(80) 



This is the desired expression for the amplitude of fluctuations at the comoving scale k2x which re-enters 
the horizon when the temperature is T2x- 

This result depends on three unknown parameters - the temperature at the end of kination Tk^end; 
and the parameters A and A' in the potential. The COBE experiment provides us with a constraint on the 
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amplitude {8p/ p{k) ~ 5 x 10^^) and the spectral index of density perturbations (0.7 < n < 1.3). The extra 
constraints we impose are those required by our consideration of electroweak baryogenesis: The phase of 
kination must persist well past the electroweak phase to have an important effect on the expansion at 
that scale. For example, for T^^ ^nd = ^ns, and T2X = leV, we find A = \/0.33 = 0.57 for A' > \/6, and 
A = \/(m = 0.33 for A' = These lead to the tilt in the power spectrum n 0.7 for A' > V^, and 
n ~ 0.9 for A' = \/5, which are consistent with the constraint from COBE. Following the discussion in 
section 4 we know that in order to create the observed BAU at a second order or cross-over electroweak 
phase transition we need to have very close to 1/a^ scaling in kination, i.e. A' > \/6- This requirement 
therefore leads in this model to a prediction of the spectral index n ~ 0.7. Using ( [79D we can also compute 
^i,end ~ -f^i,end) ^i,end) <;^i,end) ctc. In particular, for A' > \/6 wc have Ti^end = 6 X lO^GcV, and for 
A' = \/5 we have T^^ end = 2 x lO^^GeV (independent of A in inflation). 

What we have illustrated with this analysis is the observational adequacy (and even potential predic- 
tivity) of a model of this type. The 'prediction' we derived here is of course particular to a model we have 
invoked in its specific form in an ad hoc way. It would be of interest to study models which are derived 
in detail from a well motivated particle physics model. We will limit ourselves here to one qualitative 
comment on the sort of model which motivated our choice (see ) in which the parameter A varies slowly 
(logarithmically) as a function of It is not difficult to see that the constraint on the spectral index may 
be much weaker when A interpolates between our limiting values: We were constrained to increase A as 
A' increased (to give 1/a^ scaling) to make the amplitude of perturbations sufficiently large. The effect of 
having an interpolating scaling between that in inflation and scaling will be to increase -f^i^end (at 
fixed end)> which also increases the density perturbations, allowing a spectral index closer to one. 

5.4 Inflation-Kination in a Power-Law Potential 

In this section we consider another one fleld model in which an inflation-kination-radiation domination 
cosmology can be realized. Again we assume the mechanism of reheating through particle creation in 
inflation discussed in section 5.2. The potential we study is simply the non-renormalizable power-law 
potential 

''w-T^^K]^)" '''' 

where n > 4 is an integer, which is taken to be even for stability reasons. As discussed in section 5.1 this 
potential has an oscillatory solution in which the energy density in the field scales as a-6"/("+2), i.e. faster 
than radiation for n> A. It also has, as we will discuss below, "slow-roll" infiationary solutions which will 
precede such an oscillatory phase for appropriate initial conditions, just as in "chaotic" inflation in a (j)^ 
potential. As mentioned in the introduction potentials of this type have been studied in the context of 
inflation motivated by F and D flat directions of supergravity theories (see, for example, ^^). Lower 
order perturbative terms are forbidden by a discrete symmetry imposed on the superpotential. 

There are several important differences between the exponential we have considered in the previous 
section and this potential. There the potential was made up of two pieces, one with an inflationary 
attractor the other with a 'kinationary' attractor. Here we also have solutions of the two types in different 
regions of the potential, but the cross-over from one to the other is dynamically determined rather than 
an independent input (i.e. specifled by (jji^end)- Therefore once A„ and n are specified, the potential and 
^k.cnd and 6p/p{k) are completely determined. 

The second difference between the two potentials is that in the power-law potential kination is associated 
with an oscillatory mode, which can decay non-perturbatively via parametric resonance [^]. Only if such 
decay occurs after the transition to radiation domination, is the scenario we have envisaged possible. If it 
occurs a little earlier, but still sufficiently close to nucleosynthesis that the reheat temperature resultant 
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from the decay of the field is below the electroweak temperature, there will be some minor effect on the 
predictions of electroweak cosmology. We will not consider this marginal case and simply require the 
stability of the oscillatory mode until after the transition to radiation domination, which we require below 
the electroweak scale in order to have an effect on electroweak cosmology. Later in this section we will 
investigate in more detail the consequences of the resonant inflaton decay. 

In analogy to our treatment of the exponential potential in section 5.3, we now determine how the 
potential (|8^) is constrained by the requirement that (p generates cosmological perturbations of the required 
magnitude for structure formation, and that kination - driven by the oscillatory mode - ends in radiation 
domination (by the mechanism of ||l|) before the temperature T^s at which nucleosynthesis occurs. We will 
see that these two requirements cannot be simultaneously satisfied by a suitable choice of the two variables 
in the potential A„ and n. 

To determine the amplitude of the cosmological perturbations we follow exactly the analysis of the 
previous section, taking the perturbations to be given by (0). Furthermore, we will make use of ( [75|) and 
(|7^ ) to determine Ha in inflation and kination. 

The equations of motion for the homogeneous mode can be written as 



+ 3H4> + Xn 

2 



3 







1 n 



(82) 
(83) 



setting Mp = 1, i.e. with the rescaling 



6 H 
, tMp , —. 

Mp Mp 



In the standard "slow roll" approximation we take 



(84) 



2 n 



(85) 



and find from (82) that 



3 

Putting this expression back in (^) it is easy to show that the "slow roll" condition is 

b 



(87) 



Furthermore during the slow roll (p3) gives 



H din a 



n 



and hence 



a = exp , 

2n 



(89) 



where we chose for the integration constant ao = exp[—(j)Q/2n]. For completeness, we also write the solution 

12 



(n — 4)^nA^ 
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n-4 2 

t "-4 , 



(90) 
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where we chose t = such that </> ^ oo as t ^ 0. 

We can now write the desired expression in inflation: 



Ha 



— ] (p 2 exp — — (inflation) 



3n 



2n 



(91) 



which allows us to compute the first term in ([7^). 



Next we compute Ha during kination. It is convenient for this analysis to change the field and time 
variables to rescaled variables r and if given by 
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(93) 
(94) 
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The approximation of a sharp transition from inflation to kination consists in ignoring the explicit time 
dependence of Po{t), which is equivalent to 



n + 2 1 dip n + 2 

^ ~T' - 

b <p dr b 



(95) 



where ~ A„(/9q~^ is the average frequency squared of ip (see ( pTj ) below), and ipQ = (f)Q is the inflaton 
amplitude at the beginning of kination. Within this approximation, d'^a/adr'^ = (n — 4)'H^/(n + 2), and 
hence the term in the square brackets of (|9^ ) vanishes. It is this feature of the damping term in these 



variables which made their choice appropriate. Thus all of the time dependence in (|93D and (|9j) drops out 
and the equations can be easily integrated. The first integral of (p3) leads to 



Po = V{ipo) 



1 / dip\ 



2\dT 



n 



(96) 



which is just the energy conservation law for 93. 

The oscillatory solution for ip can be then expressed in terms of an elliptic integral with the frequency 
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dx 



a/1 — x^' 



where r„ is the oscillation period. Note that ui^ decreases exponentially with n, when the initial amplitude 
(/>o < 1. In the limit of a large n, Cn approaches unity. Finally, the solution to ( |94|) is 
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(where we used Hq = Hq) and 
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2(n-l) 



(kination) . 



(99) 



Hq(Iq \ CI 

The results for radiation and matter era are given by setting n = 4 and n = 2 in (^) respectively, so that 
(1?^) can be recast as 



^ix - <Pf, end = 2n Hn 



2(n-l) l' 
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The ratio of the scale factors at the beginning and the end of kination is given by ( [70| ) in section 5.2. 
Using s = 6n/{n + 2) we have (with Mp = 1) 
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and recalling that Ta = const, we have 
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so that (100) becomes 
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Taking this expression with 
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specifies the amplitude of density perturbations in the model implicitly in terms of the parameters in the 
potential and n. Comparing this to the requirement of COBE provides the first constraint on the model. 
The second constraint is the requirement that kination ends before nucleosynthesis, i.e. at a temperature 
^k, end > ^ns- After some algebra we obtain the simple relation 
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On figure 4 we show a plot of (5 = 6p/p{kix) and of as a function of n with T^^end = ^ns and 
T2X = leV. It is clear that everywhere the amplitude is too large (by several orders of magnitude) to 
satisfy the constraint from COBE. Further, to increase T^^end at fixed n we require a larger A^, which 
results in a larger 6. That kination end by nucleosynthesis thus forces the energy density in radiation at 
the end of inflation to be sufficiently large, which forces A„ to be so large that the density perturbations 
produced are too large, for any n. A single field model of this type is therefore ruled out. 

This analysis neglects the possible decay of the oscillating mode by parametric resonance into either 
its own fluctuations, or other fields that it couples to. If such decay occurs when the energy density 
of the oscillating mode has red-shifted to be sub-dominant relative to the radiation, then such decay is 
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Figure 4. 5 and A,^ vs n with T,^ end="'" 
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irrelevant and the model is simply ruled out by nucleosynthesis constraints (and the requirement that 
density perturbations not be too large). If, on the other hand, it decays when the energy density of the 
oscillating mode is still dominant over the radiation, the model may be viable if the decay products can 
thermalize with ordinary matter. In this case, however, there will always be large production of entropy 
(with corresponding dilution of the baryon to entropy ratio) and hence the model is not of much interest in 
the context of the question of how an increased expansion rate at the electroweak could lead to production 
of the BAU at that scale in scenarios when it is usually assumed to be impossible. 

It is nevertheless of interest to study the non-perturbative decay of this mode (i) to see whether such 
models are really ruled out by the observations above, and (ii) because this decay channel is relevant to a 
scenario where another field plays the role of inflaton, which we will discuss briefly below. 

We wish to compare the resonant decay time rjecay of the inflaton-kinaton to the time at the end of 
kination t; end- As we show in the appendix, the field can decay either into its own fluctuations, or to 
other fields that it couples to. In both cases the decay time Tdecay can be approximated by 

Tdecay ~ 77— ""^^ > (106) 

where no — 1/2 is the initial occupation number, and ngcatt is the 'late time' occupation number. The 
resonance is usually terminated by the back reaction effects from the created particles |45], and ngcatt can 
be estimated to be ngcatt ~ 1/Aefr = l/(-^n</'"~^nd) inflaton-kinaton decay into its own fluctuations, 

and nscatt ~ I/5 when (j) couples to a scalar field where g is the coupling constant of the interaction 
term gC,'^(fP' /2. In the appendix we also show that fi = ~ 0.16/n for the decay into its own fluctuations, 
and for the decay into another field /i ~ 0.1 — 0..2 when q > 1, and ^ <C 1 when g <C 1. 
In order to get an expression for r^^ ^nd, we make use of (H), (H), and (ITOll ) to get 

1 / 9 n 1 \^ , , 

rk,end=^- -3--^ . (107) 

l^n \ ep An (p-^ end / 
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Setting Tk _ end = Tdecay results in the following constraint on the instability coefficient 
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such that when the inflaton decays only into its own fluctuations, nscatt/''^o ~ ^ / {'^n4\~^nd) •> fJ.^^^^"^ — 
0.049, significantly larger than any fin for n > 4. This implies Tdecay — (/^n°'^***'^//^n)Tk, end, typically greater 
than Tns , where we took ^ end — ^T^s • When (p couples to another scalar field via a quartic term of form 
g(p'^('^/2, then riscatt / no ^ 1/g ~ l/(4gAefr), where q = gcfP'/Au'^ {cf. (114) in the appendix), and we see 



that for q>l, l/g < l/XcS, and hence /^^^'^ti- < 0.049. Recall that in this case /i is typically of order 0.1. 
This means that when q > 1, the inflaton decays somewhat before r^, end via parametric resonance. If, on 
the other hand, q <^ I, then ^ 1, implying a late inflaton decay, Tdecay > ^ns- 

In summary, the oscillating mode in the power-law potential decay via parametric resonance before 
nucleosynthesis if it couples to another scalar field with q > 0.1, or equivalently g > 10"^''. When g > 10"^" 
it is not immediately obvious whether the resonance shift slows down the decay or not, and, although the 
discussion in the appendix suggests that it does not, further analysis is required to establish this definitively. 
In any case we can conclude that this single power-law potential driving inflation with reheating of the 
type we have discussed (as in [^), is therefore only ruled out as a viable cosmological model for g < 10~^^. 

5.5 Two Field Models 

Finally we consider briefly models with ordinary re-heating (through the decay of the inflaton) . In this case 
the field which supports the kinetic energy dominated mode cannot also be the inflaton, but is a second 
field which comes to be important after inflation and ordinary re-heating. Again we consider the two cases 
of an exponential with its rolling mode and the power-law potential with its oscillating mode. 
Case 1: Inflaton -|- Exponential Potential 

As discussed in section 5.1 a simple exponential V{(f)) = Voe~^^^^'^'' with A such that it supports, when 
dominant over radiation, a mode scaling faster than radiation, will not come to dominate over radiation 
irrespective of the initial conditions on the field. (If the initial condition gives a scaling slower than 
radiation, it will bring the system to the attractor ( [6^ ) in which the scalar field contributes at most an 
amount comparable to the radiation.) Therefore, just as in the one-field case, a potential is required which 
is only asymptotically this simple exponential. If the field lies initially in a part of the potential which 
is flatter - flat enough to support an inflationary type solution - a period of inflation will occur once the 
radiation cools so that its energy density is comparable to that in the scalar field. The initial conditions 
and details of the potential will determine what the final ratio of the energy in the scalar field and radiation 
energy is when the scalar field enters its asymptotic kinetic mode. If this second period of inflation occurs 
at intermediate energy scales (after 'full' inflation at the GUT scale, say) and is of a small number of 
e-foldings, the ratio will be such that the kinetic energy domination may end before nucleosynthesis. 



Such a short period of inflation at an intermediate scale occurs in so-called 'thermal' inflation |50|. 
A scalar field is trapped in a false minimum by its coupling to the plasma and comes to dominate for a 
short period until the inflation it drives cools the plasma and allows it to roll away. In the present context 
all that is required is that the field, rather than rolling into an oscillating mode and decaying, roll into a 
potential which is asymptotically exponential. 

Another way in which such a transient period of inflation which cools the radiation and leaves a kinetic 
mode dominant could occur is by special initial conditions in certain potentials, e.g. if the field (p with 
potential V{(p) = Voe~^'^ I^^p sits initially close to = 0, a period of infiation will occur when becomes 
comparable to V^, the duration of which will depend on how close to (/> = the field is initially. Without 
significant fine-tuning there will be a few e-foldings of inflation followed by a period of kination. 
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And lastly, we mention a variant of hybrid inflation. Recall that in hybrid inflation, one field {(p) is 
held at the false vacuum minimum by a large expectation value of a second field {ip), and hence it drives 
inflation. When ip becomes sufficiently small, (j) rolls down to its true minimum. The roll in a steep 
potential, e.g. gcjP'ilP' /2 + V^e"'^"'^, leads to kination. Since the shape of the ijj potential determines the 
amplitude of density perturbations, we have more freedom to tune parameters of the model than in the 
one field case. In particular there is no need for variation of A. 

Case 2: Infiaton + Power-Law Potential 

The various examples just given can be carried over in an obvious way to the case of a power-law 
potential. The difference is parallel to that in the one-field case: If the field is initially sub-dominant 
relative to the radiation, the oscillatory mode about the minimum could end up being dominant depending 
on the initial conditions. If the field lies initially aX, (p > nM-p there will be a period of inflation which 
brings the field to dominance over the radiation. For a small number of e- foldings the radiation produced 
by ordinary re-heating (by decay of the infiaton) may be dominant over any radiation produced by particle 
production as in the mechanism we discussed in the one-field models. The constraints which we derived 
in the one field model in section 5.4, and which we found could not be satisfied, are circumvented simply 
because the initial radiation density is not specified by the potential, and the relation (|105| ) no longer holds. 
For the model to work we also require that the field decay via parametric resonance occur after the mode 
has become sub-dominant relative to the radiation, i.e. after nucleosynthesis, which will translate into 
some upper bound on the couplings of the field. The precise bound would have to be derived in analogy 
to the treatment given above for any particular model (which will specify end)- 

6 Testability: Connections to Other Observables 

We have shown in the preceeding section that there are simple post-inflationary cosmologies quite different 
from the standard scenario in which radiation domination begins straight after inflation. Instead a phase 
interpolates between the two in which a kinetic energy dominated mode, red-shifting more rapidly than 
radiation, drives the expansion of the Universe. As we discussed in the earlier sections of the paper, if 
such a phase continues until after the electroweak scale, the implications for electroweak cosmology would 
potentially be extremely important. To conclude we wish to discuss some of the broader implications 
of the observations we have made. In particular, we began this paper with the usual motivation for 
the consideration of electroweak baryogenesis: It promises to follow nucleosynthesis in making firm and 
observable predictions about the cosmological remnants from an epoch at which temperatures are such 
that we can have experimental knowledge of the relevant physics. It promises to be a testable theory. 
What is left of this testability now that we have effectively turned one crucial parameter, which is usually 
assumed to be known, into an unknown? 

In contrast to nucleosynthesis there is in this case only one "observable" - the baryon to entropy ratio 
- produced by a calculation. Our analysis shows that, at least in certain particle physics theories, it will 
be possible to "fit" the observed asymmetry by an appropriate expansion rate. Does making H a variable 
make the theory intrinsically untestable? The answer is negative for two reasons. Firstly, it is an extremely 
non-trivial requirement that one can produce the observed BAU in any given electroweak model, even if 
the expansion rate is a variable. In a first order phase transition, for example, the requirement of various 
parameters - most importantly on CP violating parameters - are typically extremely strong, independently 
of the expansion rate (without the sphaleron bound). As we have seen in this paper, it is conceivable that 
it could turn out that the scalar sector indicates an analytic cross-over or a weakly first order transition 
and CP violation sufficiently large that the BAU could be produced if the expansion rate is greater by 
about five orders of magnitude than usually assumed. Would we then take it this to tell us that this is the 
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case or that we are unlikely to be able to draw a definite conclusion as to whether the BAU was created 
at this scale? This brings us to the second answer to the question: The theory is truly testable only if 
we can find other observables which depend on pre-nucleosynthesis cosmology. If we do indeed find that 
the BAU can be generated with a different cosmology, that would provide a major incentive to pursue this 
possibility. 

One possibility is exactly the relic densities of dark matter particles discussed in ^ and [Q]. The 
discovery of a candidate dark-matter particle would allow one to determine the expansion rate at its time 
of decoupling from the requirement that it be the cold dark matter in the Universe. For example, from 
figure 2 in [Q] we see that the relic density of a Majorana neutrino changes by several orders of magnitude as 
the expansion rate at its decoupling does. If this indicated an expansion rate different from the standard 
value and consistent with that required at the electroweak scale for generation of the observed BAU, 
one would have compelling evidence that cosmology is indeed different. Another possible way of probing 
pre-nucleosynthesis cosmology is with magnetic fields, which in certain models are produced at or before 
the electroweak scale, or at the QCD phase transition. This seems a more remote possibility for a firm 
constraint in that the connection to observed fields is itself very indirect. However, it is one worth bearing 
in mind. For example, in the mechanism discussed in |31] in which fields are generated by an instability 
related to the abelian anomaly, the expansion rate enters in determining when perturbative processes come 
into equilibrium. This depends on the expansion rate, and for a significantly different expansion rate the 
results would be different. 

Further there is also the possibility of probing cosmology at the electroweak scale indirectly by its 
connections to other epochs. A good example of this is in fact the scalar field cosmology we have discussed, 
in particular the exponential potential. In this case the same coherent mode which dominated in kination 
can in fact play an important role again at later times. We noted the existence of an attractor solution 
with energy densities given as in ( |68[ ) for the exponential in the presence of a component of matter or 
radiation. How soon this will be established after the end of kination depends essentially how much the 
ratio of kinetic to potential energy at the end of kination differs from its value in the attractor (|6^), and 



this will vary depending on the model. In |11| the case is treated in which this transient period between 
the two attractors is assumed to end well before the beginning of matter domination, and details of the 
observable consequences on structure formation in a flat CDM dominated universe are studied; in |12] 
the case of entry into the attractor well into the matter era at a red-shift z ~ 70 is treated. With the 
forthcoming satellite experiments which will measure the properties of the microwave sky, such models will 
become testable in detail. 



7 Appendix: Resonant inflaton decay in the potential 



In this appendix we study the decay of the infiaton via parametric resonance. First we address the decay 
into its own fiuctuations, and then we discuss how it decays into other fields. We start with writing the 
evolution equation for small fiuctuations around the infiaton zero momentum mode: (f) ^ (j) + 6(f) in (|82|). 
After a Fourier transform and setting M 
equation 



Mp = 1 as in (84), we get the following linearized mode 



+ 3H 



dt 



+ 
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Rescaling to the new variables as in ( |9^ ) and assuming pure kination, i.e. that the field amplitude is small 
), we obtain 



4(n-4) 



5v9g = 0. (110) 



Notice the scale dependence next to k"^ which means that, even though the zero-mode equations are time 
independent, the mode equations are not. Assuming adiabatic variation of a/ao, Eq. ( |11C| ) becomes the 
famous Hill equation 



^= — — + g = -r^ , 111) 

where r' = u;„t and / = 2(v?(r')/v?o)''"^ - 1 is defined so that max|/| = 1, (/) = 0, /(r' + tt) = /(r'). 
The general solution of Hill's equation is of the form e^^'^ P{'^')-, where P(r' + tt) = P{t'), and it is often 
given as the stability chart. The unstable regions are specified by the curves of constant n in the {q, A} 
plane, and the stable regions are bounded by /i = 0. The instability chart is important since the field 
decays exponentially into the modes with > 0, preferably so to the ones with large /x. The special case 
of the Hill equation when n = 4 - the Lame equation - is extensively studied in the literature on inflaton 



decay |53, 54], [55|. The instability chart exhibits unstable regions which branch off from A = n? at q = 0. 
For q <^ A one is in the narrow resonance regime, since the bands are narrow and /i <C 1. The chart 
is symmetric under q —q. On the other hand, for 1 < 2q < A the resonance bands become broad 
and fi "large." Typically, when 2g ~ ^, peaks at ~ {2tt)^^. In this case the field decays very fast, 
characteristically in a few dozens oscillations. 

Notice that in general for a given n, g„ = (n — l)(/9Q~^/(4a;^) is specified. Consequently, to get a 
rough estimate of the decay time, it suffices to plot the one dimensional slice g = g„ of the chart. As 
the field decays, q stays constant, unless the backreaction of the created particles is large enough to 
change ujn- Numerical simulations []52| , |48[| show that for the X^ip^^ potential the backreaction from created 
particles grows to about 6m'^ = 3X{6ip'^) ~ A49J0/4, changing the effective frequency a;^ — > u;^ + drn^, and 
consequently reducing q to about half of its original value and A ^ A + dm? /{uJ^ + 5m?'). The growth 
of dm? is terminated by narrowing the resonance as a consequence of the backreaction on A and g, and 
intensifying scatterings of the resonant particles off the zero mode, as a consequence of increasing resonant 
amplitudes. By then a significant portion of the field has decayed. We will assume that a similar scenario 
holds for a generic 0" case. This is plausible since, as we will see below, the instability charts are quite 
similar. 

We have evaluated numerically [^] the instability charts for some of the models. The results for 
n = 4, 8, 16, 32 are plotted in figures 5. The corresponding initial values for g„ are 1.045, 3.84, 14.29, 54.64. 
Note that in all cases to a very good approximation the first instability band terminates at A = 2q, so 
that the field decays into the second instability band. The values for ^ are /i4 = 0.0425 0, //8 = 0.023, 
file = 0.011, = 0.0056, so that /x„ ~ Afi^/n ~ 0.16/n. For the following analysis the details of the chart 
are not that important. It is sufficient to keep in mind that for n larger, /i^ decreases. 

As discussed in section 5, an inflaton that decays into its own fluctuations and does not couple to other 
fields leads to disastrous consequences for nucleosynthesis. Indeed, since the inflaton decay products scale 
as radiation and eventually dominate the energy density of the Universe, and also decouple from the rest of 



-"Note that this value is a bit higher from 1x4 — 0.0359, the value quoted in 
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Figure 5. iJ-^ vs A 
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matter, they behave effectively as many additional massless degrees of freedom, leading to a very different 
expansion rate than predicted by nucleosynthesis. If, on the other hand, the inflaton predominately decays 
into another scalar field, which consequently thermalizes, producing thus standard radiation and matter 
particles, nucleosynthesis may be unaffected by late inflaton-kinaton decay, as long as the decay occurs 
comfortably before nucleosynthesis. 

Before we start discussing the inflaton decay time, we outline the physics of the inflaton decay into 
another scalar field C- We assume a standard quartic coupling g to C (that itself couples to Standard Model 
particles) of the form ^(^^0^/2 such that the linearized mode equations of motion are 



+ 5HCk + 



a 



0. 



(112) 



With the rescahng and C = C(ao/a)'^/("+2), this becomes 



+ 



4(n-4) 
n + 2 



6(n-4) 
n+2 



9^ 



(113) 



and can be recast as 



A. 



C 



dT'2 ^ 

4(n-4) 
n + 2 



6(n-4) 
n+2 gf^ 

4u;2 



(114) 



^ In this case conformal rescaling might seem more appropriate since it would get rid of all dependence on a. Nevertheless, 
we stick to the rescaling in (p2|), in order to be able to make direct comparison of decay times. 
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where r' = a;„r and = 2{lp{t') / Lpo)"^ — 1 is defined so that max|/^| = 1, (/^) = 0, /^(t' + vr) = /({t'). 



As above in (111), in adiabatic hmit, this reduces to Hill's equation. There are however two differences: 
first, can assume a wide range of values depending on g, and, second, is a (growing) function of a. 
The corresponding Mathieu equation for n = 2 is studied in great detail in the literature, and shows that 
for > 1 the field decays with an average value ~ 0.1 [ p7| . For the conformal case with n = 4 a similar 
value for fi is obtained. Here we will assume that, for any n > 4, ^ ~ 0.1 as well. 

Now we present an estimate of the decay time Tdecay For a moment we assume that the resonance 
shift does not drastically affect particle production. Later on we comment on the plausibility of this 
assumption. The typical initial mode amplitudes are such that the corresponding initial "occupation 
numbers" tx uJk'Pj^ip_j^ (where ujk is the energy of the mode k) are of order jt,™'*'^' ~ no — 1/2. Since 
the resonant mode amplitudes grow as 6ipk oc exp fii^ojnT, one can estimate the field decay time as follows. 
The field decays when the energy in fiuctuations become comparable to the energy in the zero momentum 
mode, i.e. when the occupation numbers ~ no exp 2^i^ujnT become of order ~ 1/Aeff (Aeff = An</>o~^), 
implying that the decay time can be approximated by 

Tdecay ~ TT" • (US) 

2a;„/ifc no 

This same equation applies for the field (f) decaying into other scalar fields. The only difference is that 



the maximum occupation number is in this case ngcatt ~ 1/5- As a caveat to ( |115| ), the authors of |48] 
showed that one should expect longer decay times if the self-coupling of the second field (C) is large, i.e. 

^ g. In this paper we do not dwell on these complications, and assume the couplings such that the 
simple estimate ( |115| ) is valid. 

Finally we comment on how the time dependence of A in ( p.llD and ( |114| ) can affect the decay time 
( |115| ). We first discuss the decay into a second scalar field. The (comoving) resonant momentum is specified 
by 5 A ~ yfq, which in (114) gives 



4 



2 y/g4'Q^n / Oo \ "+2 



'ci.-^^^j . (lie) 

This agrees with the well known result that for n = 4 the resonance is static. On the other hand, for 
n = 2 (Mathieu case) the resonant momenta grow rather fast as the Universe expands and, for q large, 



adiabatic approximation breaks down, leading to "stochastic resonance" [57|. However, it turns out that 
the instability exponent is rather robust and maintains the value /i ~ 0.1. What happens when n > 4? 
In this case the resonant momentum decreases with time and again for g > 1 we expect breakdown of 
adiabatic approximation. Just like in the n = 2 case we expect /i to be robust and be of order /x ~ 0.1. 
This should not in any case be considered as proof, but conjecture. 



In the case when the field decays into its own fluctuations, 5 A ~ ^Jq gives (c/. ( Ill )) 



4(n-4) 

kL ^ \^{n-l)Xn <^'?u^n (^) (117) 

which again leads to a shift in the resonant momentum. Unfortunately, the conformal case (n = 4), in 
which the resonance is static, is the only case studied in the literature, so we cannot make any analogy 
as we did in the former case. Since in this case the resonance is rather narrow, the resonant momentum 
redshift may significantly slow down the decay. An implication would be that the effective /i decreases, 
leading to somewhat less stringent bounds on A„ and n than indicated in (|108|). 
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